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CSJ I In the present paper we describe the set of form factors for hadronic r decays based on 

Resonance Chiral Theory. The technical implementation of the form factors in FORTRAN 
ijT^ ■ code is also explained. It is shown how it can be installed into TAUOLA Monte Carlo 

^ . program. Then it is rather easy to implement into software environments of not only Belle 

• ! and BaBar collaborations but also for FORTRAN and C++ applications of LHC. The 

Q I description of the current for each r decay mode is complemented with technical numerical 

CN I tests. The set is ready for fits, parameters to be used in fits are explained. Arrangements to 

work with the experimental data not requiring unfolding are prepared. Hadronic currents, 
ready for confrontation with the r decay data, but not yet ready for the general use, cover 
/\^ '. more than 88 % of hadronic r decay width. 
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1 Introduction 

Measurements of r lepton, because of its long lifetime, large mass and parity sensitive 
couplings, lead to broad physics interest. From the perspective of high-energy experiments 
such as at LHC, knowledge of r lepton properties offers an important ingredient of new 
physics signatures. From the perspective of lower energies, r lepton decays constitute an 
excellent laboratory for hadronic interactions. In itself, the r lepton decays constitute an 
excellent laboratory for studies of hadronic interactions at the energy scale of about 1 GeV, 
where neither perturbative QCD methods nor chiral Lagrangians are expected to work to 
a good precision [U |2l |3l IH |5] . At present, hundreds of milions of r decays are amassed by 
both Belle and BaBar experiments. It is of utmost importance to represent such data in a 
form as useful for general applications as possible. 

Most of these data samples are not yet analyzed. For example, in Ref. [6j only 10 % 
of the collected sample, which means 5.4 M events for r^ — )■ z/^vr^vr^, was used. Future 
samples at the Belle II or Frascati Super B facilities will be even larger [71 [8] . That means 
that already now the statistical error for the collected samples of r — > SttUj. is of the order 
of 0.03 %. For r -^ Kmru^ it is about 0.1 % and for r -^ KKnu^ at the level of 0.2 %. To 
exploit such valuable data sets, theoretical predictions need to be properly prepared. As 
typically several milions of events per channel are collected, that means that the statistical 
error can reach ~0.03 %. To match it, parametrizations of hadronic currents resulting from 
theoretical models must be controlled to technical precision better than 0.03 % in Monte 
Carlo, combining theoretical aspects and full detector response. Only then, one can be sure 
that the comparison of the data with theoretical predictions exploits in full the statistical 
impact of the data, and one can concentrate on systematic effects both for theory and 
experiment. One should stress that the above technical precision is required not only for 
signal distribution, but for background as well. Sophisticated techniques allowing proper 
comparisons of the data and models are also needed. 

A review of the status of available tools for such studies of r decays can be found 
in Ref. [H]. It was concluded in Ref. [9] that the appropriate choice of hadronic current 
parametrization was the most essential missing step to perform. It was also found that, 
for the decays involving more than two pseudoscalars in the final state, the appropriate 
use of hadronic currents in fits is important too. At present, standards of precision are at 
2 % level. This is a factor of 100 less than what is required. For many r decay modes even 
this 2 % precision level is far to be reached [TO] . 

The original version of TAUOLA [TT] uses the results of Ref. [121 [13] and their extensions 
to other decay channel^j. In that model each three-pseudoscalar current is constructed as a 
weigthed sum of products of Breit-Wigner functions [3 [121 [IHl [I6]. This approach was con- 
tested in Ref. [1^ where it was demonstrated that the corresponding hadronic form factors, 
which were written to reproduce the leading-order (LO) xPT result [18j, fail to reproduce 
the next-to-leading-order one (NLO) [T^ [20] . The corresponding parametrization based 



^With time, due to pressure from the experimental community, many other parametrizations were 
introduced, but not in a systematic way. Some of those found its way to TAUOLA later, Ref. [T3], and are 
used as a starting reference point for our present project as seen from the computing side. 



on Breit-Wigner functions was not able to reproduce CLEO r^ — > {KKTT)^i'r decays data 
|21] . This resulted in the CLEO collaboration reshaping the model by the introduction 
of two ad-hoc parameters that spoilt the QCD normalization of the Wess-Zumino part. 
This shows that, although the approach of weighted products of Breit-Wigner functions 
was sufficient and very successful twenty years ago now, with the massively increased ex- 
perimental data samples, it is pressing to upgrade. As an alternative, an approach based 
on the Resonance Chiral Theory [221 123] was proposed. Its application to hadronic tau 
decays is supposed to be consistent and theoretically well founded (see Sect. [7] for the 
related discussion). However, its results have to be confronted with the experimental data 
before actual improvement will be confirmed. The hadronic currents for the two and three 
pseudoscalar final states that we consider here have been calculated in the framework of 
RxT [m 1211 [251 |26l [27] and have been prepared for TAUOLA. 

Section |2] is devoted to a general presentation of the form hadronic currents must fulfil 
to be installed into TAUOLA generator ^llj. In each subsection analytic forms of currents 
calculated within Resonance Chiral Theory are given channel by channel. In Section [3] 
energy-dependent widths as used in the parametrization of intermediate resonances are 
presented. Section Hj is dedicated to technical tests of the channel r — t- SirUr- For channels 
involving kaons, only overall benchmark distributions are collected. The details of technical 
tests are left to the project Web page |28j. For all decay channels numerical results, 
which are of more physical interest, are collected in Section [5l Within it, the three-meson 
channels, which have been worked out in more depth, are first presented and then, the two 
meson channels are discussed. In both cases technical aspects are worked out to precision 
better than 0.1 %. 

The organization of the hadronic currents and how they can be integrated into the 
TAUOLA library is explained in Section El and in Appendix [Bl Section [7] is prepared for 
a reader who is oriented towards the theoretical details of calculation and estimation of 
theoretical uncertainties of the approach. It provides arguments necessary for discussion 
of the range of parameters allowed for fits. The summary in Section [8] closes the paper. 
Further technical appendices are also given. Appendix |X] lists analytic functions used in the 
parametrization of hadronic currents. Appendix O provides numerical values of the model 
parameters used all over the paper. It is explained which parameters and in which range 
can be modified without breaking assumptions of the model and where in the code they 
are defined. Appendix [D] collects branching ratios of the newly prepared r decay channels 
as calculated by Monte Carlo simulation. It contains explicit references to the definition 
of the hadronic current in each decay channel, which are spread all over the paper. In the 
future, however, these definitions can be replaced by the new references. 

The implementation of final state interactions (FSI) in the two-meson r decay modes 
is discussed in Appendix [El An improvement of the latter and the scalar form factor in 
r — )■ KtcUt decays will be addressed in a forthcoming publication |29] . 

Finally let us stress that our paper aims at explaining how this new set of hadronic 
currents can be installed in TAUOLA, independently of whether it is a standalone version, 
part of Belle/BaBar software or a different configuration, which is another purpose of 
Appendix [Bl 



2 Hadronic current for two and three hadrons 

Before discussing in detail the implementation of the currents into the program and re- 
sulting distributions, let us first collect here all necessary formulas. In general we will 
follow conventions for normalizations as used in Ref . [11] . We will not recall here relations 
between hadronic distributions and decay product distributions though. They are rather 
simple and we assume that the reader is familiar with the necessary parts of Ref. |11] . 
Let us recall the matrix element Ai for the r decay into hadronic state X and a neutrino: 
t{P) -)■ Xiy^{N). It reads M = ^u{N)Y{l-^r^)u{P)J^. All dynamics of hadronic inter- 
actions is encapsulated in a current J^, which is a function of hadronic r decay products 
only. 

Contrary to TAUOLA, as documented in Refs. [H] or [T3], now hadronic currents for all 
two-pseudoscalar final states are defined in separate routines and the constraint that the 
scalar form factor must be set to is removed. 

For r decay channels with two mesons [/ii(pi) and /i2(P2)]; the hadronic current reads 



J^ = N 



P1-P2- — (Pi +P2) ) F^'is) + ^{p,+p,YfP{s 



(1) 



where s = {P1+P2Y and A12 = ml — ml. The formulas for vector, F^(s), and pseudoscalar, 
F^{s), form factors depend on the particular decay channeo and are given, respectively, 
for the tt'tt", {KTr)~ and K~K^ decay modes in subsection 12.41 and the following oneCl- 
SU{3) symmetry relates all four normalization factors by the appropriate Clebsch-Gordan 
coefficient: 

72' V2 2 ^ ^ 

For the final state of three pseudoscalars, with momenta pi, p2 and ps, Lorentz invari- 
ance determines the decomposition of the hadronic current to be 

J^^ =N {Ti^[c^[p^-p^yF^ + C2{p:,~PirF2 + C:,{p^-p2rF:,] 

+c,q^F, - —i^cse^^^KPsPs^s}, (3) 

where T^y = g^y — g^gi^/g^ denotes the transverse projector, and q^ = (pi + P2 + Ps)'^ is 
the momentum of the hadronic system. The decay products are ordered and their four- 
momenta are denoted, respectively, as pi, p2 and p^. Here and afterward in the paper F 
stands for the pion decay constant in the chiral limit. 

Functions Fi (hadronic form factors) depend in general on three independent invariant 
masses that can be constructed from the three meson four-vectors. We chose q^ = {pi + 



^The vector form factor of both two pions and two kaons is expected to be fixed at zero momentuni 
transfer by gauge invariance in the SU{2) symmetry hmit ^U\: F^(0) ~ 1, see Section [^^ 

•^Two-meson r decays involving an 77 meson, t ^ rj^ 'P^v-r, P — tt,K have a negligible branching 
fraction l3n[32l. 



P2 + Pd)"^ and two invariant masses Si = {p2 + Ps)^, S2 = {pi + Ps)^ built from pairs 
of momenta. Then S3 = (^1+^2)^ can be calculated from the other three invariants, 
S3 = q"^ — si — S2 + ml + m\ + m|, and Fi written explicitly with its dependencies reads as 
Fj(g^, Si, S2). This form of the hadronic current is the most general one and constrained 
only by Lorentz invariance. For modes with an even number of kaons the normalization 
factor reads as A^ = cos6'cabibbo/-^5 otherwise A^ = sin6'cabibbo/-^- 

We leave the F4 contribution in the basis, even though it is of the order ~ rn^/q^ [T7t[25l 
[15] for the three-pseudoscalar channels that we consider in this work. It plays a role in the 
low q^ region for the three-pion modes. We will neglect the corresponding contribution for 
the modes with kaons, i.e., C4 = 0. Among the three hadronic form factors which correspond 
to the axial- vector part of the hadronic tensor, (Fi, F2, F3), only two are independent. We 
will keep the definition of Fi, F2 or F3, exactly as shown in Eq. ([3]), which is the form used 
in TAUOLA since the beginning. However, linear combinations constructed from only two of 
these functions are in principle equally good. The decay channel dependent constants Cj 
are given in Table [H 



Decay mode 

(Pl, P2, Pz) 


Cl 


C2 


C3 


C4 


C5 


TX^Tl^Tl^ 


1 


- 1 





1 





n%'7i- 


1 


- 1 





1 





K-7V-K+ 


1 


- 1 








1 


K%-K^ 


1 


- 1 








1 


ir-TrOi^o 





1 


- 1 





- 1 



Table 1: Coefficients for formula ([3]) in the isospin symmetry limit. Note, that in Ref. [TT] 
different conventions were used and coefficients were affecting normalization too. 



The theoretical assumptions behind the hadronic currents that we use are discussed in 
Section [3 In the model the results for all hadronic currents, with the exception of two- 
pion and two-kaon modes, are calculated in the isospin limit, therefore the corresponding 
hadronic form factors depend only on the average pion [m^ = (m^Q + 2 ■ r?7,^+)/3] and kaon 
[rriK = {rnxo + m^+)/2] masses, we relax the assumption latei|j. For the three-pseudoscalar 
modes every hadronic form factor consists of 3 parts: a chiral contribution (direct decay, 
without production of any intermediate resonance), one-resonance and double- resonance 
mediated processes. The results for the hadronic form factors are taken from Refs. [211 125] . 

The channels we will present (together with the trivial decay - from the Monte Carlo 
point of view - into ttVt and Kvj.) represent more than 88 % of the hadronic width of r 
[3T]. The dominant missing channels r — )■ Ti^TT'ir'TT^Vr and r — )■ tt'ti^ti^tt^Vt-i which are 
together about 9.7 % of the hadronic width of r, are more difficult to control theoreticalljo. 



"^At first step we will take such assumption for our phase space generator as well. Later in the paper we 
will nonetheless return to proper masses, distinct for charged and neutral pseudoscalars. We will evaluate 
the numerical consequences, see Tabled 

^ Several older options developed for these channels are provided for user convenience, but they will 



Also attempts to describe r — ^ ii^ii^Ti^Ti^i'r are relatively recent. They are technically 
compatible with our solution for TAUOLA currents and can be used simultaneously. It is 
documented in Ref. [33] . 

Let us now describe hadronic currents for each particular channel. 



2.1 7T 7T TT^Ur and 71 7T 7T 



0_0_- 



l^r 



Hadronic form factors for the three-pion modes have been calculated assuming the isospin 
symmetr}o, as a consequence, Tn^± = m„o . The code for the current is given in 
new-currents/RChL-currents/f 3pi_rcht . f . 



The independent set of hadronic form factors Fi. 



5 is chosen as: Fi, F2 and 



Fi (-F3 = then). The vector form factor vanishes for the three-pion modes due to the 
G-parity conservation [16ll33]: F5 = 0. It is convenient to present the functions as 



(F^ + Ft + Fn-R 



377 



1,2,4 



(4) 



where F^ is the chiral contribution, F^ is the one resonance contribution and F^^ is the 
double-resonance part. The R^'^ constant equals -1 for 7r°7r*^7r^ and 1 for 7i^7i^tt~^. 

For the convention defined by Eq. ([3]), the form factors Fi can be obtained from Ref. 
with the replacements: 



Fi{Q\s,t)^Fi{q^,Si,S2)/F, ^ 

where F was defined after Eq. ([3]). The form factors read: 

2^2 



1,2,4, 



F^{q\s,,S2) 
Fl{q',SuS2) 



(5) 



(6) 



2G 



V 



V 



V2FvGx 
3F2 

2g2 _ 2si 



3si 



Sl 



Friq',s,,S2) 



M2 
3F2 g2 



s, - Ml - iMpVp{si 
+ 



p ---pi p^ 
■S3 , S3 



Sl 



iMj:,{si 



•'P'- p\ 



S2 



Mj 



q^ 



Ml-iMATAiq") 



- iMpTp{s2) 
- (A' + A") 



3 Sl 



+H 



£1 ^4 



2g2 + Sl - S3 



Sl 



M'p 



iMpTpisi] 



+ H 



£2 ^ 

9 ' 9 
qz qz 



Sl 

S3-S1 



M2 



iMpVp^Si) 



S2 



Mj 



iMpTp{s2) 



not be documented here. Please see the README files stored in directory new-currents/other-currents 
for details. 

^The inclusion of the complete first-order corrections to SU{2) symmetry is beyond our present scope. 
In particular, electromagnetic corrections arising at this order are neglected. We restrict ourselves, for 
the moment, to the ones given by the mass splittings between members of the same S't/(2)-multiplet. 
These enter the kinematical factors and phase-space integrals. The model-independent electromagnetic 
corrections can be handled with PHOTOS [Ml El], while the structure-dependent corrections have been 
computed only for the one- and two-meson r decay modes (with only pions and kaons) in Refs. [SHJ [SZl 
|38l 139]. Its implementation in generation with PHOTOS will follow work of Refs. [40l |4T1 142] . 



where 

H{x,y) = ~\y + \'x + \\ (7) 

and 

X' 



2V2FaGv ' 

4Ao = A' + A". (8) 

Bose symmetry implies that the form factors Fi and F2 are related F2{q'^,S2,si) = 
Fi{q'^, 81,82) (the minus sign that comes from the definition of the hadronic current, Eq. 
([3]), is included in C2). 

The pseudoscalar form factor, F4 = F^ + F^, carries the contribution from both the 
direct vertex and the one-resonance mechanism of production: 

2 - _ 2V2 ml[3{83 - ml) - q\l + 2kR^^)] 



F^iq ,81,82) 



3 2g2(g2 — m^) 

Fliq"^, 81,82) = 3^2 ^ [Q2(g^^S2,^Sl) + a2(g^gl,^S2)] , (s) 

where k = 1 for r^ — ;■ tt^tt^tt+z/^, k = 1/2 for r^ — )■ Tr^vr'^Tr^z/^ and 



/2 ^ BGy gj m^ ^3-^2 .,„^ 

The pseudoscalar form factor F4 is proportional to m1/q^ [2l], thus it is suppressed with 
respect to Fi and F2. However, the pseudoscalar contribution can affect the q^ spectrum 
near the thresholqll. 

Besides the pion decay constant F, the results for the form factors Fi depend on some 
coupling constants of the model: Fy (we impose Gy = F^/Fy), Fa and the masses of 
the nonets of vector and axial-vector resonances [My and Ma) in the chiral and large- A'c- 
limits. We follow Refs. p3| 125] and replace the masses used in the resonance Lagrangian 
with the masses of the corresponding physical states: My — )■ Mp anqj Ma — )■ M^^ . 

To include the p' meson we follow Ref. p3], its Eq. (32). We insert Eq. fITT]) of the 
combined p and p' propagators into our Eqs. ([6]) and (fTOl) 



1 1 



M2 _ ^2 _ ^MTJq^) 1 + /3„ 



^ Pp 



Mj - g2 _ zM,r,(g2) Ml, - g2 _ iMp,Tp,{q 



(11) 



^Numerical results with and without i^4 are presented in Section |4?2] 

^ In footnote ^^of Appendix [0.2 1 we explain; the parameter Ma of the short-distance QCD constraints 
should not be identified with Ma^ ■ Our choice Ma -^ Ma^ is not well founded, but can be easily changed 
at the time of fits to experimental data. 



Impaclo of the p' meson on the dV /dq^ spectrum can be seen from Fig. 3 of Ref. [24] . 

In the file new-currents/RChL-currents/f 3pi_rcht .f the form factors Fi, F2, F4 
given by Eqs. ©-(iQl) and with substitution (fTTI) are coded. For completeness, let us 
remember that form factors F3 and -F5 are equal to zero. 

The only eventual isospin breaking will result from mT^± ^ ^.,^0 used in the phase space 
generator embedded in tauola.f. 

2.2 K-TT-K+Ur and K^n-R^Pr 

Again isospin symmetry is assumed, therefore m.,^± = m^a and m^^o = m^±. The code for 
the currents is given in new-currents/RChL-currents/fkkpi .f . 

We will neglect the contribution from the pseudoscalar form factor F4 as it is propor- 
tional again to rn^/q^ [25]. We present the result for the non-zero form factors Fi in the 
same way as before: 



FJ" + FT + Ff 



1,2,5. 



(12) 



Taking into account the convention for the current, Eq. ([3]), and the values of the q 
coefficients in Table [U the result for the form factor can be obtained from Ref. [25] with the 
replacements: Fi{Q'^,s,t) -)■ ^2(5^, S2, Si)/F, F2{Q^,s,t) -)■ Fi{q^,S2,si)/F. Therefore, 
the form factor Fi reads 



FHq' 
F-{q' 

Friq' 


S2,Si) = 
S2,Si) = 

S2,Si) = 


V ^ 

3 ' 

V2 FvGv 
6 F2 

2FaGv 

3 F2 Ml- 


[ M2 -S2- iMpTp 

q' 
- g2 _ iMATAiq") 
A««(g2, 


(^2) Ml.-s,- 
"5^^(g^si,S3,S2,m 


(1 
iMk*Vk'{si) _ 

9 9 9 \ 


[ M2 -S2- iM, 

si -iMK*TK*{,si) 


rp(s2) 



where the functions A", B^\, A^^ and B^^ are defined in Appendix |X] and S3 is calculated 
from q^,Si^S2 and masses. 



^A discussion on the implementation of the second (third) resonance nonet to the model can be found 
in Section [T] 



The form factor F2 is given by 






Fril^ ^2, s,) 



px 



(14) 



V2FyGv 

6 



F2 



A«(g^S2,S3,m^,m2,mf^) 



+ 



2 FaGv 

3 F2 



MI.-s^-iMk*Tk*{si) 



M2 - S2 - ^M,^,(S2) 



Ml 



g2-^MAr^(g2 



M2 



S2 - iMpTp{s2 



+ 



i?^^(g^S2,S3,Sl,mLm 



2 m|f) 



M|. 



si -iMK*TK*{sx) 
The vector form factor, F5, arises from the chiral anomaly [HI HS] and the non- 



anomalous odd- intrinsic-parity amplitude |16]. It is obtained from Ref. [25] with the 
replacement Fs{Q'^,s,t) — )■ —F5{q'^,S2,Si)/{47r'^F^). It reads 



F^{q'^,S2,Si 
Fliq^S2,Si 

+ COS^ 6'y 



(15) 



Wti^Gv 
1 -y/2tan9v 



^R/ 2 2 2 2\ I ■ 2 n 

C [q ,S2,mK,mK,mJ I sm Oi 
C"^(g2,si,m|^,m2,m|^) 



l + y2cot^y 

M2 - S2 - iiv^r. 



2Fv 



M2 - S2 - zM^r 

^"(g',S2,Sl) 



Gy M2 - g2 _ ,M,r,(g 




M|. -Si-zM^,r;^.(si) 



^KK( 



, S2,Sl, 



-16y27r2FyG' 



y Ury 



C-(g^32,m^) 



sin^ ^i 



M2 - g2 _ zMpr,(g2) 

l + V^COt^y 2/, 

' COS &i 



C™(g2,si,m|^) 



Mi, -Si-2M^,rx*(si) 



1 - v^tan6'v 



M2 - S2 - zM^r, 



Mi-S2- iMsT, 



where C", -D^ and C^^ are defined in Appendix |Al 

For the widths of the narrow resonances u and the PDG [31] values are taken and the 
constant width approximation is followed. The parameter Oy defines the mass eigenstates 
a;(782) and 0(1020) [25j and is the mixing angle between the octet and singlet vector states 
Wg and Uq. In our numerical calculation we take the ideal mixing \0y = tan~^(l/\/2)]. In 
this limit, the contribution of the 0(1020) meson in Eq. flT5|) vanishea^l. This will be 
changed in future with fits to the data. 

The file new-currents/RChL-currents/fkkpi .f contains the form factors Fi, F2 and 
F5. The form of Eqs. ( TT3|) -(TT5l) is used. As one can see, the contribution from the excited 
states, e.g. p', p", is not included in the KKii case, contrary to the three-pion one. The 
only eventual isospin breaking assumed will result from m^± 7^ m^o and mK± 7^ rrixo used 
in the phase space generator embedded in tauola.f . 



^^ On the other hand, one should keep in mind that the decay channel r 
BaBar [47] and r" -^ (j)K~Vr by Belle [l8] . 



Att v-r was observed by 



10 



2.3 K-TT^K^Ur 

The hadronic current of the K'tt^K^v^ decay mode is again obtained with isospin sym- 
metry; m^± = m^o and irixo = mx± . The code for the currents is given in the file 
new-currents/RChL-currents/f kkOpiO . f . 

If the momenta of the pseudoscalars are attributed as in Tabled! that is K~{pi)n^{p2)K^{p3) 
then it is convenient to choose the independent set of hadronic form factors for the axial- 
vector part as -^2(0'^, -52, S3) and F^^q"^, 82,83). We will neglect the contribution due to the 
pseudoscalar form factor F4 as it is again proportional to the square of pion mass over 
q"^ [25]. Taking into account the constants q of Table [Hand Eq. ([3]), the result of Ref. [25] 
needs the replacements 

F2{Q',8,t) -^ F3(g^S2,S3)/F, Fl(Q^s,t)^F2(g^S2,S3)/F, 



F3{Q\8,t) 



-)■ 



As before: 
and 

F^iq^,82,83) 
Fl{q^,82,83) 



-F,{q\82,83)/{An^F'' 



F]" + Fr + Fr 



IFvG 



v^v 



6 F2 



B^{82,8i,m 



Ki 



m: 



(16) 
2,3,5 (17) 

(18) 
A^(g^S2,Si,m^,m2,m^) 



Ml,-83-iMK*rK*{83] 



+ 2 



M2 



82 - iMpTp(82 



+ 



A«(g^Si,S2,m2,m^,m 



K) 



Ml. 



Friq',82,83) 



V2FaGi 



81 -iMK-'TK*{8i) 

5^^(g2, S2, si, S3, m^, ml, m\) 



3 F2 Ml-g2_,M^r^(g2) 



Mh 



83 -iMK*TK*(83) 



+ 



^ A''''{q'^,82,8i,mj^,ml,mj^) ^ A^"(g^ si, 82, ml, m\, m\) 



M2 - S2 - lMpVp{82) ' Ml, - Si - iMk'Vk'{8i 

83 is calculated from q^, si, S2 and masses. The contributions to F3 read 

S2,S3) = 0, 



F^"' 



(19) 



Fl{q^,82,83) 



IFvG 



yu-y 



6 F2 



A''{q'^,83,8i,ml,ml,ml}_ ^ 5«(s3, Si,m|^, 



m 



K) 



Ml 



83 - iMk*Tk*{83) ' " M2 - S2 - lMpVp{82) 

A"(g2^ Si, S3, rn\, m\, ml) ' 

Ml,-8,-tMK^TK'{8,) \ ' 
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FriQ\s2,s, 



V2FaG^ 



3 F2 Ml-q^-tMATAiq^) 






+ 2 



5««(g2, S3, Si, S2, m^, mf^, m^) A«« (g^^ si, S3, m|^, m|^, ml) 



M2 - S2 - ^Mprp(s2) M|., - si - «M^.r^.(si) 

The form factor F5 driven by the vector current is given by the sum of 
F^{q',s,,ss) = 0, (20) 



F^{q\s2,s,) 



8V2tt^ Gi 



C^(g2, S3, m^, ml, m\) C'iq^, Si, m\, ml, m\) 



2Fv E^{s3,si) 



Friq',S2,ss) 



Gy Ml - g2 _ ^M,r,(g2 



-IQti'^FvG 



V^V 



Ml 



C^^{q^,S2„m 



K) 



MI,-S3-iMk*Tk*{s3) 



C^^{q^,Si,m 



K) 



MI-.-s^-iMk*Tk*{si) 



where the new function E^ , as the previous ones, is defined in Appendix \K\ 

The file new-currents/RChL-currents/fkkOpiO.f contains the form factors F2, F3 
and F5. Eqs. flT8|) - fl20l) are used. As one can see, the contributions from the excited states, 
e.g. p', K*' are not included in the form factors of the r — )■ KKt^Vj. decay channels. This 
shall be an obvious future improvement resulting from the confrontation with the data. 

The assumption that the eventual isospin breaking will result only from m^± 7^ m^^o 
and m,K± 7^ triKO as used in the phase space generator embedded in tauola.f is taken. 



2.4 tt^ttV^, TT'^R-Ur, TT^i^V^ and R-R'^Vr 

The two-pseudoscalar final states are simpler. They can be presented together in one 
subsection. The code for the hadronic currents of the tt^tt^z/^, K~K^Vt- and {Ktt)"!/^ 
modes is given, respectively, in files new-currents/RChL-currents/frho_pi .f , 
new-currents/RChL-currents/fkOk.f and new-currents/RChL-currents/fkpipl .f . 

In the general case there are both vector and scalar form factors. In the isospin sym- 
metry hmit, m^± = m^o, m,K± = ^n^o, for both two-pion and two-kaon modes the scalar 
form factor vanishes and the corresponding channel is described by the vector form factor 
onl}o. Also for the Ktc mode we restrict ourselves at first to the vector form factor only, 
the scalar form factor will be properly included in Ref. [29] . 



^^Thc scalar form factor appears only at ncxt-to-leading order in SU(2) breaking and it can be safely 
neglected (see Ref. [35] for details). 
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For all three channels we use a parametrization for the vector form factor which is 
developed starting from the lowest-lying resonance contribution: 



fU-'^) 



-:VMD, 



s) exp 






loop 



967r2F2 



ReApQ{s) 



(21) 



where F^^^ is the contribution from the lightest vector resonance that can be exchanged 
in the process, and the exponentiation resums FSI effects (see Appendix |E] for a related 
discussion). The function Apq[s) is a loop function for two pseudoscalars with masses nip 
and niQ, and is presented in Appendix |X1 N^J^ is a constant dictated by chiral symmetry 
for the different decay channels: 



Nr^° = 1 Ar,^"^° 

loop ' loop 



r\1 loop 



nZ 



3 

4' 



(22) 



In the case of the two-pion mode, the theoretical calculation performed in the framework 
of RxT [19] gives the following result for the form factor: 



F^As) 



Ml 



Ml-s-tMpTp{s] 



■ cxp< 



967r2F2 



RcAt^-^o^s) H — ReAK~Ko{s) 



, (23) 



if only p resonance is taken into account. The contribution of the excited resonances (both 
p' and p") modify the form factor of formula ( l23l) . Following Ref. |19] it takes the form 



FLis) 



M2 + s(7e*'?^i +(5e^'^2) 



—s 



M^-s 


- iMpTp{s) 


- exp 


M}- 


s - iMpiTpi 


w" 



967r2F2 



ReA^-„o{s) H — ReAx-Ko{s) 



-sVp, 



^Ml,al{Ml,] 



Ml„ -s-iMpnTp„{s) 
where the phase-space factor 

ap{q') : 



exp< 



-sVpn 



7,Ml„al{Ml„] 



ReAJs) 



ReAJs) 



(24) 



1 — Am]:>/q^ 



(25) 



is used. The function Ap is the same loop function Apq defined in fH5|) . but in the limit of 
equal masses. It is given in Appendix Rl In our file new-currents/RChL-currents/f rho_pi . f 
the form factor of formula fl2^ is used. 

The two-kaon vector form factor is written following Ref. [2S] as 



M2 
p 



KK\ 



M^-s- iMpTpis] 



exp 



{ 



967r2F2 



ReA^^^o{s) + -ReA^-Rois) 



(26) 
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One can see the expression for the kaon vector form factors coincides with the pion form 
factor from Eq. (l23l) . 

The excited resonances have not been taken into account in Eq. ( |26l) following Ref. |26j . 
However, their implementation along the lines of Eq. f[Ml) is simple. The file 
new-currents/RChL-currents/fkOk.f contains botlo forms of the two-kaon form factor 
dUD and (ESD. 

For the Kn mode we applied the result of Ref. [27] Eq.(5), which reada^^l: 

pv /„^ ^ ( Ml, + s^K. sjK. A . . 

^"^ ^ \Mj.,-s-iMK*TK4s) Ml„-s-iMK*'TK*,{s)) ^^ 



—s 



exp< 



1287r2F2 



ReAxTris) + ReAx^is) 



Note that due to the FSI effects, the form factor F^^(O) 7^ 1. 

The m^± 7^ m^o and mj^± 7^ 171^0 used in the phase space Monte Carlo generator, 
and discussed numerically later in the paper, will be the only isospin breaking assumed. 
However, we plan to include electromagnetic corrections to r — t- -K^ir^Ur decays in the 
future (seel^. 

3 Energy-dependent widths of resonances 

In this section we collect the formulas to calculate the energy- dependent width of the 
resonances p, p', p" , K*, K*' and ai. They were used in the previous section as ingredients 
for the construction of hadronic currents. From the technical side their calculation requires 
integration of the appropriate matrix elements over the phase space. In this way, for 
example, unitarity constraints are taken into account |2H |5T1 152] . 

The energy-dependent width of p(770) resonance, calculated in the SU{2) limit (m^± = 
m^rO, rriK^ = rriKO ), is given |52] as 



r iq') = ^^ 



al{qy{q' - Aral) + I^WMq' " 4m^) 



f28) 



The phase-space factor a-,, is defined in Eq. ( l25ll . the 27r and 2K loops are included. This 
form of the p width is used in modes of three pseudoscalars. 



^^By default, our program runs with the K^K~ vector form factor of Eq. (pS)) . However, changing the 
value of the parameter FFKKVEC = to FFKKVEC = 1 in value_parameter .f allows to run the code with 
the form factor of Eq. (f24|) . Numerical effects due to the inclusion of the excited resonances are given in 
Section EH 

^■^By default, our program runs with the Kt: vector form factor of Eq. ([27]). However, changing the 
value of the parameter FFKPIVEC = 1 to FFKPIVEC = the code will run with the form factor given in 
Eqs. (17), (18) of Ref. JSOj. For discussion, see Section [Ol 
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In the two-pion and two-kaon modes SU(2) breaking effects were taken into account in 
the p off-shell width 



rM' 



967rF2 



e{q' - thr^^)X'/' (l, ^, ^) + le{q' - thrKK)\"' (l, 



2 2 



q^ J I \ q'- q'- 

where A(x, y, z) = {x — y — zY — %z, t/ir^7r = (^-,7+ + m^o)^ and thrxK = ("^x+ + ^/^o)^. 
At this stage the widths of the p'(1465) and p"(1700) mesons are modeled as decays to 
two pions, 

with r,. = r,. (Af2) . 

The energy-dependent width of the i^*(892) resonance is given, in the SU{2) limit 
(r?7,7r± = m^o, mK± = m^o), in Ref. [53] • It is related to rp(g^) by chiral symmetry. It 
reads: 



,(29) 



TM^l' 



1287rF2 



,3/. (^1 !^_ ^^ ,(,. _ j^,^j ^ ,3/. (^i_ !^_ ^^ ,(,. _ ,,^^^.j 



,(31) 



with thrxn = {rnx + Tn.„Y and thrxr] = ij^K + ^r?)^- 

Formula fl3Tl) is used for the KKtc modes whereas for the Kir modes we use the 
following result (Eq. (4) from Ref. [27] ) 



TK^g' 



r 



K* 



j2 A3/2 (1, < , f^ e{q' - thvK^) + A3/2 (1, !^, <) 9{q' - thr^,) 



Mh 



A3/2(l 



-I 



-' mJ7' m^. 



+ A3/2(1 



-' mJ7' m^. 



From Eq. (4) of Ref. [27] we have 



r^. = r^*(M|.; 



647rF4 



A=^/2 (^1 



771 



ii" 



niz, 



+ 



\"^ Ti, 



m 



ft: 



r2 ' 



mz 



Ml. 



(32) 



(33) 



however, we prefer to write down Eq. (1321) in terms of the width (F^*), see discussion after 
Eq. (17) in Ref. [27]. The width of i^*'(1410) is modelled as a decay to Ktt and reads 



A3/2 f 1 



K- 



'{<f) 



m\ ml 



K* 



Mk*' A3/2 f 1, < 



_m2_\ 



■e{q^ - thrKn) 



For the energy dependence of the ai resonance widtlo we usqlj 
^aAq') = 2Tl{qy {q' - Qml) 

+ 2Vl\q')e {q^ - (m. + 2mK?) + Ff (g^)^ {q^ - (m. + 2mKf) 



(34) 



(35) 



^^Calculation of oi width from Fi was already used in Refs. [5l fT2l[54] . 

^^There is an additional factor 1/F^ here and in Eq. (|38|) compared with the definition, e.g. in Refs. [24l 
[2S]. This is related to the normalization of our form factors Fi, F2 and F3. See explanation prior to the 
Eqs. di), dni) and ([T8|). 
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where 

"' ^^ ' 192(27r)3F2F2M„, V (P 

J dsdt {VfF, + V^F2 + ViF^Y^'' {{Vi^Fi + V2^F2 + Vs^FsT'"')* (36) 

stands for the contribution from the individual three-pion and (two kaons - one pion) 
absorptive cuts. Here 

V:' = c,T^''ip,-Pk)., t^j^k=l,2,3, (37) 

with the coefficients q appearing in Table 1. 

In summary, r^^(g2) is the contribution of the tt^tt^tt'^ and tt'^tt^tt^ cuts, F^ (g^) of the 
K^-K^K^ and K^-k^K^ cuts, and finally the K^tc^K^ contribution gives rise to the term 
r^ (g2). The form factors Fi are presented in Sections 12. 1^ 12.21 and I2.3[ The symmetry 
factor is defined as S" = l/n\, where n denotes the number of identical particles in the final 
state. 

For reference, we include the formula for the spectral function, the g^-spectrum for the 
processes r — )• 3 pseudoscalars Ut of this worlO 



dq 
where 



dq^ 128(27r)5M.F2 V g2 ') J «^«^ 



1 /^. , o 9 



WsA + -{l + 2^]iWA + WB) 



,{3l 



Wa = -iVl'F, + Vi'F2 + V,^Fs){Vi^Fi + V2^F2 + V3^Fs)\ (39) 

WsA = q'IFaI'. 

The following phase space integration limits have to be used 

t+{s) 



dsdt = ds dt, (40) 

where 

t±(^) = ^{(g^-m^)^-[A^/^(g^^,m^)TA^/^(m^,^,m^,^,.)p}. (41) 

The necessary functions are located in files new-currents/RChL-currents/f unct_rpt . f 
and new-currents/RChL-currents/wid_al_f it .f . 



^^Our testing programs feature a calculation of spectral functions, including those of the formula (p9|) . 
However, we will not elaborate on this point here, even though it is important for future data analysis [9], 
where results of Ref. [13] are proposed to be used. We expect that a Monte Carlo sample will be used 
instead of semianalytical Eq. (p8|) . 

16 



4 Benchmark calculations for three-pion mode 

Since Ref. [H] has been published, numerical tests of TAUOLA Monte Carlo functioning 
have not been repeated in a systematic way, despite the technical precision requirements 
are much higher now and reach sub-per mil level. Prior to physics oriented comparisons 
between analytical (numerical integration) and Monte Carlo calculations, we need to re- 
visit numerical stability of the generator and of multiple numerical integration, used in 
semianalytical calculations accompanying generation and its tests as well. 

We will use the decay channel r — )■ tt'tt^tt^Uj. to demonstrate our tests of Monte 
Carlo. For other channels, technical tests will be skipped from document at ioio even 
though new issues absent in the r — )■ tt^tt^tt^z/^ case can appear. A good example is 
numerical stability at phase space edges of presamplers for relatively narrow resonances 
such as K* in r — > KirUr decay. In this case, for long runs square root of a negative number 
may have appeared because of rounding errors; an appropriate correction to the code was 
introduced. Such long runs were never performed in the past for this channel. Our tests 
were indeed long. In some cases, for variables counting crude events we have even run over 
the allowed maximum (~ 2 ■ 10^) for FORTRAN INTEGER type. 

To avoid problems with multidimensional integration of the ai-meson propagator which 
is rapidly-changing as a function of its arguments, we first tabulated the r^^ (g^lof Eq. (^^. 
The code for tabulation is located in new-currents/RChL-currents/tabler |^. Then we 
use linear interpolation to get the value of the ai width at required g^. To integrate over 
s and t variables the Gauss integration method has been used. The produced distribution 
has been checked to be numerically stablqlj. 

4.1 Technical test 

Before we can go to the presentation of simulation results, where physical currents are used, 
let us start with the simplified cases. We will begin with the calculation for r — )■ 7r~7r~7r"'"i/^ 
where Fi = F = 0.0924 GeV, other form factors are set to zero {F2 = 0, F4 = 0) and 
m^± = m^o = 0.13804 GeV. It is an important starting point, it helps to adjust conventions 
of normalization constants in TAUOLA Monte Carlo and analytical calculations. Phase space 
integration is free from singularities resulting from the matrix elements. Nonetheless, 
corresponding presamplers can be verified. Numerical integration is rather quick as there 
is no need of invoking time consuming functions. That is why it is important to perform 
this check with a precision higher than for later tests. 



^^Directory new-currents/RCiiL-currents/cross-check is devoted to such tests. The README file 
explains technical details. 

^^Technical details of the calculations, which are quite independent from parts of the code loaded with 
TAUOLA, are explained in README files of this directory and its subdirectories. 

^^Thc source code of the integration routine gauss has been taken from CERN program library [55] . 
Tests are provided in new-currents/RChL-currents/cross-checks . Cross-check by linear interpolation 
of the q^ distribution from the neighbouring points demonstrates that the fluctuations due to numerical 
problems of integration are absent, the results are continuous, whereas the result produced with the 
integration method VEGAS [SS] had a tendency to fluctuate. 
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The total rate we obtained from a Monte Carlo run of 6-10^ events was (2.7414±0.01%)- 
10~^^ GeV. For seniianalytical numerical integration we obtain (2.7410±0.02%)-10^^^ GeV 
(with 2- 10~^ precision tag). A difference of 0.015% was found. This technical test performs 
better than could be achieved at time of work for reference [H]. Statistical samples are 
larger by 3 orders of magnitude than what could have been used at that time. 

For completeness, let us provide a plot of J^ generated from Monte Carlo divided by 
the semianalytical (numerical integration was used) result for this spectrum. Reasonable 
agreement is found, see Fig. [H 

Triple Gaussian integration is used for the analytical calculation and double Gaussian 
integration enters as well into the current calculation to be used in matrix elements of 
Monte Carlo generation. That is why pretabulation for the ai width, raj(g^), is convenient 
as it speeds generation enormously. This represents another technical feature being tested 
by normalization study and figures like Fig. [H 

Technical tests, as the one we discuss now, belong to the group of comparison booklets 
collected in the Web page 



4.2 Test with semirealistic parameters 

Let us now introduce the physical content of the current, keeping at first m.,^± = 771^^0 = 
0.13804 GeV and dropping out statistical factor | of two identical vr's. Agreement between 
Monte Carlo and semianalytical calculation should be ~0.01% again. 

We take Fi as given by Eq. (Q, but all other currents are set to 0. Our numerical 
results for the rate from numerical integration, (1.8721 ± 0.02%) ■ 10~^^ GeV, and from 
Monte Carlo generation, (1.8722 ±0.01%) ■ 10"^^ GeV (run with 6 • 10^ events), agree well. 
The difference is only 0.005% thus compatible with statistical error of generated sample. 
Comparison of differential distributions, analogue to Fig. [H is available from the Web 
page [28]. For the next step we choose the Fi and F2 form factors according to Eq. ([6]), 
F4 = 0. Numerical results are: (4.2015 ± 0.02%) • 10^^^ GeV for semianalytical calculation 
and (4.2023 ± 0.01%) ■ lO'^^ GeV for Monte Carlo generation. The difference is 0.03%, as 
expected. The figures are again available from Ref. [2B]- 

Finally, we consider the result for the total width predicted by Resonance Chiral Theory, 
given by the Fi, F2 and F4 contributions in Eqs. (|6]) and (|9]). The F4 form factor does not 
affect the value of the width at our precision and in this case it is (4.2025 ± 0.01%) ■ 10~^^ 
GeV (run with 6 ■ 10^ events). 

The above comparisons check also that the differential distribution dV /dq^ in r — ?■ 37ri^^ 
is numerically stable. In particular that it is not affected by numerical problems due to 
double Gaussian integration. We have also checked, that if the function value has been 
obtained from interpolation of neighboring g^'s, the difference with the value calculated 
directly was appropriately small as should be expecteco- 



'^^ The appropriate program for a test is available in the directory 

new-current s/RChL-currents/cr OSS-check/ check_analyticity/check^nalyt_3pi 



18 



2.0 



1.5 



1.0 



0.5 - 



0.0 



TTT — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — r 



Ratio: MC/semianalytical, Fi = F 




J_LJ I I I I I I L. 



_] I I I I I I I L. 



_] I I I I I I I L. 



LJJ I I I I I I L. 



1.0 



2.0 



3.0 



4.0 



'TT^U^ 



Figure 1: An example of a figure included in tests collected in project web page 
The ratio of the histogram and analytical formula is shown for the J^ in r — t- tt^tt^: 
decay mode. Only Fi = F is non-zero. Normalization is not adjusted. A statistical 
sample of 6- 10^ events was used and semirealistic initialization as explained in this section. 
Agreement within statistical errors is found. Fluctuations at the ends of the spectra are 
due to substantially less populated bins of that region. These tests represent a technical 
test not only of Monte Carlo generation but also for semianalytical numerical integration. 
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Channel 


Width, [GeV] 


PDG 


Equal masses 


Phase space 
with masses 



K-TT-K+ 



(5.778 ± 0.35%) ■ 10 
(9.72 ± 3.5% ■ 

(1.9 ±5% ; 
(3.60 ± 10% ; 

(2.11 ±0.8% 
(2.10 ±1.2% 

(3.17 ±4% ] 

(3.9 ±24% ] 
(3.60 ± 12.6% 



10-13 




10-15 




10-14 




10-15 




10-13 




10-13 




10-15 




10-15 




■ 10-15 





(5.2283 ± 
(8.3981 ±0 
(1.6798 ±0 
(2.0864 ± 
(2.1013 ±0 
(2.1013 ±0 
(3.7379 ± 
(3.7385 ± 
(2.7367 ±0 



,005%) 
,005%) 
,006%) 
,007%) 
,016%) 
,016%) 
,024%) 
,024%) 
,025%) 



10-13 

10-15 
10-14 

10-15 

10-13 
10-13 
10-15 
10-15 
10-15 



(5.2441 ± 
(8.5810 ±0 
(1.6512 ±0 
(2.0864 ± 
(2.0800 ± 
(2.1256 ±0 
(3.8460 ± 
(3.5917 ±0 
(2.7711 ±0 



,005% 
,005% 
,006% 
,007% 
,017% 
,017% 
,024% 
,024% 
,024% 



10-13 

10-15 
10-14 

10-15 
10-13 

10-13 
10-15 
10-15 

10 
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Table 2: The r decay partial widths. For each channel, the PDG value [3T] is compared 
with numerical results of Monte Carlo integration of our currents. The third column 
includes results with isospin averaged masses, whereas for the last column physical masses 
were used. Comparison of the last two columns enumerates the numerical effect of physical 
masses, breaking the assumption of isospin symmetry in a potentially uncontrolled way. 
Further results for individual decay channels are given in Subsections of Section [5l 

5 Numerical results for two and three-pseudoscalar 
channels 

In the previous section we have presented examples of technical tests. Let us now con- 
centrate on numerical results, corresponding to the most refined options of the currents 
included in our distribution tar-balo which are of physics interest. 

In the phase space generation, we will take into account the differences between neutral 
and charged pion and kaon masses, physical values will be taken. This has to be done to 
obtain proper kinematic configurations. On the other hand, this choice breaks constraints 
resulting from isospin symmetry (see l^ in a potentially uncontrolled way. That is why we 
collect numerical results from Monte Carlo calculation in the form of Table [21 where the 
partial widths from Particle Data Group compilation |3I] are compared with our results 
obtained with isospin-averaged pseudoscalar masses and with the physical ones. 
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Figure 2: The r — )• 7r~7r~7r"'"i/T- decay: comparison of distributions for TAUOLA cleo cur- 
rent [14] and for our new current. On the left-hand side, the plot of 7r~7r~7r"'" invariant 
mass is shown and on the right-hand side n~^7i~ invariant mass is given. Green histograms 
(light grey) are for the new current, red (darker grey) are for TAUOLA cleo. Distributions 



for the r — t- 7r*^7r°7r z/^- 



decay coincide with the ones for r — t- tt n n 



-^-^+, 



5.1 



TT TT TT 



Ur and TT^TT^TT Ur 



Let us now turn to numerical results for the r — t- Sttz/^ decays obtained with our currents 
of Subsection 12.11 From the semianalytical calculation of partial width we get (2.10073 ± 
0.02%) • 10-13 QgY fo^ ^ ^ ^_^_^+^^ ^g^^y^ ^^^ (2.10072 ± 0.02%) ■ lO'^^ Q^y ^^ 

T — )■ 71^71^71^ Ur, practlcally the same value. The approximation of the equal masses for vr^ 
and 7r° has been taken. In this case the Monte Carlo results are identical for tt-tt-tt'*" and 
7r°7r°7r- final stateS we have obtained (2.1013 ± 0.016%) ■ lO^^^ geV. 

For the physical, i.e., distinct m,r± and mT^o, masses we have obtained (2.0800±0.017%)- 
10-13 GeV for 7r-7r-7r+ and (2.1256±0.017%) ■ lO'^^ GeV for tt'^tt^tt- mode. The difference 
for the distributions is too small to be seen and we present plots for the tt^tt^tt- case. Only 
two example plots are given in Fig. [21 We point the reader to the web page [2S], for the 
booklet of comparisons obtained with MC-TESTER [ST]. The figures for dV /dq^ spectrum 
from Monte Carlo and analytical calculations are also available from the plots of the Web 
page [28]. 

From the technical point of view to separate generation of the two 37r sub-channels one 
has to set the BRAl = for the Tr^vr'^Tr- mode, and BRAl = 1 for the Ti^n^Ti^ , (e-g-, in 
routine INITDK, which is defined in our demonstration program 
new-currents/Installation/demo-standalone/taumain.f ). 

An attempt on comparisons of the new model distributions and experimental data is 
also given later in the paper, in Section 15. 5[ 



^^For the convenience of updates, we have prepared Appendix IdI to be modified in the versions of the 
present paper to be included in the tar-ball. It summarizes elementary installation benchmark results, 
that is the branching ratios calculated by the Monte Carlo. 

^^The chiral contribution F^{q^, si, S2), Eq. ([9]), differs for tt^tt" 
does not affect sizably the width. 



However, Fn{q^ , S2, si) 
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5.2 K-Ti-K+v^ and K^tt-R^ 



l^T 



As in the case of Svr decay modes, the Monte Carlo generated distributions are relegated 
to the project web page [2S]. In particular successful checks with the analytic function 
for dV /dq^ taken from Ref. [25] are shown there. In the following, we present figures |3] 
and m comparing the two histograms obtained with our new and cleo versions of TAUOLA 
currents. The differences are substantial. CLEO Collaboration [58| [2T] introduced an 
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Figure 3: The r — )■ K~tt~ K^Ut decay: comparison of distributions for TAUOLA cleo 
current [H] and for our new current. On the left-hand side, plot of K~7t~K^ invariant 
mass is shown and on the right-hand side K~^7t~ invariant mass is given. Green histograms 
(light grey) are for the new current, red (darker grey) are for TAUOLA cleo. 
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Figure 4: The r — )■ K^Ti~K^Vj. decay: comparison of distributions for TAUOLA cleo cur- 
rent [14] and for our new current. On the left-hand side, plot of K^gK^Ti~ invariant mass is 
shown and on the right-hand side K^ti~ invariant mass is given. Green histograms (light 
grey) are for the new current, red (darker grey) are for TAUOLA cleo. 

ad-hoc parameter violating a property stemming directly from QCD, the normalization of 
the vector form factor given by the chiral anomaly [59]. Of course, in our new current we 
are not taking into account excited resonances. In the past, two couplings could only be 
estimated because of unavailability of data to determine them from a fit. These aspects 



22 



can (and should) be improved at the time of confronting with the datacj. 

Let us now turn to the decay widths. The result for SU{2) symmetric masses from 
TAUOLA with a sample of 2 ■ 10^ events is T = (3.7379 ± 0.024%) ■ 10"^^ GeV for K-7r-K+ 
and r = (3.7385 ± 0.024%) ■ 10"^^ GeV for K%-K^. The difference for the partial width 
of the two channels is within statistical error. The analytical result is the same for both 
channels and was found to be (3.7383 ± 0.02%) ■ 10^^^ GeV. It agrees with the ones of 
Monte Carlo. 

For physical masses of the pseudoscalars, the Monte Carlo results for K~7r~K~^ and 
K%-K^ are, respectively, T = (3.8460 ±0.024%) ■ lO'^^ GeV and T = (3.5917 ± 0.024%) • 
10^^^ GeV. The effect of the mass adjustment in phase space gives an effect of the order 
of 3%. Agreement with PDG results (see Table [2]) is acceptable, but improvements from 
fits are envisaged. 



5.3 K-n^K^u 



Again, Monte Carlo generated distributions are relegated to the project web page 
In particular, successful checks with analytic function for dV /dq^ taken from Ref. [2S] are 
shown there. Figure [5] presents comparison between histograms obtained with present 
and cleo versions of TAUOLA initialization. As we can see, the differences are substantial. 
Explanations of the previous subsection for the K^tt^K^ case apply. 



xlO^ 





Figure 5: The r — )■ K'TT^K^Ur decay: comparison of distributions for TAUOLA cleo cur- 
rent [14] and for our new current. On the left-hand side, plot of K~t[^K^ invariant mass is 
shown and on the right-hand side K~tt^ invariant mass is given. Green histograms (light 
grey) are for the new current, red (darker grey) are for TAUOLA cleo. 

Let us now turn our attention to the partial width for K~7r^K'^. The SU{2) symmetric 
result from TAUOLA and a sample of 2 ■ 10^ events, (2.7367±0.025%) ■ 10"^^ GeV, agrees well 
with the analytical result (2.7370 ± 0.02%) ■ 10^^^ GeV. The effects of realistic masses in 



^•^Whenever possible, the agreement with the data should not be achieved by straightforward violation 
of the theoretical assumptions. Discrepancies may point to faulty background subtraction, or call for 
improvements or replacement of the model used in currents calculation. In practice, this may be difficult 
and require significant and simultaneous effort on both theoretical and experimental sides. That is why 
one may have to accept temporary introduction of ad-hoc factors into the currents now as well. 
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phase space are of the order of 1 % and the Monte Carlo resuh is (2.7711 ± 0.024%) ■ 10"^^ 
GeV. Agreement with the PDG result (see Table [2]) is not good. 

Analytical results for K~7t~K^ and K~tt^K^ ^ which are obtained in the SU{2) limit 
r(r- -^ K-TT-K+iy)/T{T- -^ K--n^K^u) = 3.7383/2.7370 ~ 1.366 ~ 4/3, compare well 
with the result of Ref . [25] . 

With this subsection, we complete a presentation of results for r decays into three 
pseudoscalars. Let us now turn to the case of the decays into two pseudoscalars, which are 
simpler from technical point of view. 

5.4 tt^tt^Ut-, t^^K^Ut-^ TT^K^Ur and K^K^Ut 

In this case, there is only one non-trivial invariant mass distribution, dT/ds, which can 
be constructed from the decay products. This distribution and its ratio to semianalytical 
result is given for all two pseudoscalar final states in the Web page j28]. For all two- 
pseudoscalar modes we use samples of 2 ■ 10^ event a^^l. 

The analytical result for r" -^ i^^n-n'^ equals (5.2431 ± 0.02%) ■ 10"^^ GeV, and for 
T- -^ UrK-K^ is (2.0863 ± 0.02%) ■ 10"^^ GeV. The Monte Carlo results are respectively 
(5.2441 ± 0.005%) • 10-13 Q^y ^nd (2.0864 ± 0.007%) ■ 10"^^ GeV. In both channels the 
physical values of pion and kaon masses are used. As one can see, the obtained K~K^ 
width is only ~ 58% of the PDG value. Since the mass of the p resonance is less than 
the two-kaon threshold, a significant contribution has to be expected from both p' and p" . 
To check this assumption, we used the parametrization (12^ for the two-kaon form factor 
Prk- For the moment, we use the same numerical value of the parameters 7 and 5 as in 
the two-pion case, an assumption which holds in the S'f/(3) symmetry limitcj. The result 
for the partial width of r" -^ p^K'K^ is (2.6502 ± 0.008%) ■ 10"^^ GeV. However, in the 
real world the parameters for pion and kaon modes are likely not to coincide and have 
to be fitted by the experiments. Corrections of order ~ 30% due to SUiZ) breaking are 
expecteco. The Monte Carlo result for the sum of the two channels 7i^K~Ut- and ti~K^i'^ 
is (2.5197 ± 0.008%) ■ 10"^^ GeV if the SU{2) symmetric masses are used and a sample of 
2 ■ 10^ events is generated with TAUOLA and new currents. The corresponding analytical 
result is (2.5193 ±0.02%) ■ 10"^^ GeV. The TAUOLA run with physical pion and kaon masses 
gives (2.5092 ± 0.008%) ■ 10"^^ GeV. Separate partial widths for -n^R-Ur and -k-R^u. 
channels calculated from the Monte Carlo are given o in Table [2l This can be compare^^^' 



^^By default, we include the FSI effects, the parameter FFVEC = 1. FSI can be switched off if FFVEC = 
is set in the file new-currents/RChL-currents/value4Daraineter .f . 

^^The parameters 4>i and 4>2 are subleading and their values are unsubstantial for this check. 

^^For the two-kaon mode one can include a contribution from excited resonances by setting FFKKVEC 
= 1. To run with only p(770) exchange, Eq. ^, set FFKKVEC = 0. Our defauh is FFKKVEC = and 
FFVEC=1 . 

^"^To separate submodes one has to set BRKS = for -k'^K^ or BRKS = 1 for 7r"^° in routine INITDK 
residing, e.g., in our example new-currents/Installation/demo-standalone/taumain.f . 

^^ To run the Kt: mode with the vector form factor of Eq. ([271), FFKPIVEC = 1 should be set in 
new-currents/RChL-currents/value4Daraineter.f . To use Eqs. (17) and (18) of Ref. [SO] FFKPIVEC 
has to be set to 0. For the default we take FFKPIVEC = 1. 
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with the result 2.1829 • 10^^^ GeV for the KirUr partial width based on Eqs. (17) and (18) 
of Ref. [50]. The numerical value for parameters of our model are taken from Ref. [50] , 
Table 4, second column. The differencCj between our and [50] is aboulo 15% for the Kn 
partial width. 

In order to test our improvement in the treatment of FSI lj we have run the program 
with the FFVEC = 0. It corresponds to neglecting the real part of the loop contributions in 
FX(s), F}^j^{s) and F^^(s), namely ReApgis) = 0, ReAp{s) = in Eqs. (^-(^. In this 
case, the results for the partial widths are T{t~ -)■ z/^7r"7r°) = (4.0642 ± 0.005%) ■ 10"^^ 
GeV, r(r- -^ VrK-K^) = (1.2201 ± 0.007%) ■ lO^^^ GeV (note that only the p meson 
exchange is included, namely FFKKVEC = 0), T{t- -^ v^-k^R-) = (7.4275 ±0.004%) • 10"^^ 
GeV, r(r- -^ UrTT-K^) = (1.4276 ± 0.006%) • 10-^^ GeV. 

Comparison with the data is technically simpler and further phenomenological efforts 
should be delegated to that stage of the work. In fits, one has to work with single one- 
dimensional dr/ds spectrum and current as a function of single argument s as well. We 
present comparisons of dT/ds spectra from cleo current and the present parametrization of 
TAUOLA, for n^TT^Ur, K" K^Vr in Fig. |6l and for TC^K'Ur and tt^K^Ut in Fig. [TJ respectively. 





_x10'' 
900 
800 
700 
600 
500 
400 
300 
200 
100 



0.2 0.4 0.6 0.8 1 1.2 1.4 1.6 1.8 2 



0.2 0.4 0.6 0.8 1 1.2 1.4 1.6 1.8 2 



.--,^0 



TT^'h'r and T ^ K K^Ur decays: Comparison of distributions for TAUOLA 





Figure 6: r ^ n 

cleo current [H] and for our new current. On the left-hand side, plot of tt'tt" invariant 
mass is shown and on the right-hand side K~K^ are for new current, red (darker grey) are 
for TAUOLA cleo. 



5.5 Attempt at comparison with the data 

Let us stress once again that our parametrization for all new currents is based on Resonance 
Chiral Theory and is thus self-consistent. However, only minimal attempts on adjusting 
to the data have been performed. Only one dimensional q^ distributions have been used 



^^We are thankful to Jorge Portoles for the discussion on the differences between the models. 

■^^ The reason for this difference of 15% is because we are not using Eq. (19) of Ref. 50 . The difference 
between the results of Refs. [27] and [50] for the vector form factor contribution is only ^ 4%. See Section[7] 
and Appendix |E] for a related discussion. 

^^Even a simple Breit-Wigner includes a crude description of FSI [60], where, e.g., off-shell cffcts are 
neglected in the resummation of loops. 
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Figure 7: r ^ 7r^K~i>r and r — ;■ 7r~ K^Ur decays: Comparison of distributions for TAUOLA 
cleo current [H] and for our new current. On the left-hand side, plot of tt^K~ invariant 
mass is shown and on the right-hand side vr~i^^ invariant mass is given. Green histograms 
(light grey) are for the new current, red (darker grey) are for TAUOLA cleo. 



in case of Svr [21] and KKn |2S] decays, to fit parameters such as Fy, Fa etc. For the Svr 
channel, relatively good agreement with ALEPH data is shown in Fig. 3 of Ref. [21], but 
it represents a consistency check of the input. The proper work on fits is only to start now. 
The computing and theoretical framework is ready. It is not surprising that, for example, 
agreement with the unfolded BaBar data, see Fig. |8]of Ref. [61], is not satisfactory. 
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Figure 8: Invariant mass distribution of the tt+tt^ pair in r — t- tt^tt^tt+z/^ decay. Lighter 
grey histogram is from our model, darker grey is from default parametrization of TAUOLA 
cleo. The unfolded BaBar data are taken from Ref. [61] . The plot on the left-hand side 
corresponds to the differential decay distribution, and the one on the right-hand side to 
plot ratios between Monte Carlo results and data. Courtesy of Ian Nugent. 
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6 Program organization 

The reference version of the TAUOLA hbrary that we used as a starting point for the present 
work is the TAUOLA cleo documented in Ref. [H]. The choice is not accidental. This 
version is used as a starting point for Belle and BaBar software as well. Also C++ im- 
plementation of TAUOLA [62], installed in the library of Ref. [63], uses this initialization. 
Porting into collaboration software is thus technically prepared. This means the solution 
will be convenient for many users. 

Hadronic current represents a rather compact segment of the simulation package. That 
is why only minor changes need to be introduced to TAUOLA cleo code and the makefile. 
Only several lines will be necessary to modify in the collaboration software. All the rest is 
included in tar-ball to be expanded in tauola directory. See Appendix IB] for details. 

The user will be able to switch to the new current invoking simply CALL INIRChL(l) 
prior to TAUOLA initialization, still retaining possibility to use the old ones with CALL 
INIRChL(O) as well. 

An algorithm for working with auxiliary weights to implement simultaneously several 
models of r decays is a straightforward extension. 

6.1 Weight recalculation 

Present day experimental data feature very high precision over all directions of multidi- 
mensional phase space. Nonetheless cross contamination between different channels takes 
place. This is the case, for example, if particular decay channels differ by presence or 
absence of vr^'s. Figure 6 of Ref. [6] represents such folded comparison between data and 
Monte Carlo for the r+ — )■ tt+tt'^z/^ decay channel. This result can not be used, without 
additional information on other decay channels, for fits of hadronic current. This takes 
place even for this seemingly simple case where hadronic current can be directly deciphered 
from a one- dimensional distribution. 

In general, for the precision matching of data and models it is convenient to simultane- 
ously confront several models and take into consideration all decay channels simultaneously. 
Such a solution may be a necessary technical step if fits for unfolded data are envisaged at 
the precision level better than few percents. 

To facilitate technical tasks, a weight recalculation method is prepared, following dis- 
cussions and recommendations given in Ref. [9]. It was agreed that organization of the 
programs should enable simulation following the scheme of Fig. |9l For the single generated 
Monte Carlo sample (when all detector and experimental acceptance effects are taken into 
account), one can calculate weights enumerating change of matrix element. The procedure 
can be repeated as many times as needed. In this way the optimal choice is met when 
all experimental and theoretical effects can be taken into account. With this, the fit to 
multidimensional distribution of measured data can be determined. 

We have prepared the necessary changes for TAUOLA. The following algorithm was 
checked to work (its technical details are given in Appendix IB.4P : 
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Figure 9: Weight recalculation model, as discussed in Ref. [M] . 



1. For each generated r stored in a datafile, the user program reads flavours and 4- 
vectors of r and its decay products. 

2. Appropriate kinematic transformation is performed enabling a recalculation of the 
matrix element, exactly as at the generation step. 

3. An appropriate routine of TAUOLA is chosen and the matrix element is calculated. 

4. Another instance of TAUOLA initialization can be activated and the matrix element 
recalculated using a different physics model. 

5. If one of the two models was used in the generation of the user sample, then the 
weight for model replacement can be calculated as a simple ratio of the two. 
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Usage of this method of weight recalculation does not require any changes to TAUOLA 
code except those related to installation of the presented upgrade, aiming at installation 
of our new RChL currents. Needless to say, the method can be used several times for user 
variants of RChL and interpolation of weights can be used in fits. In our example, given 
in directory new-currents/Installation-Reweight we provide a simple method to read 
generated events from file, but it can be easily adapted to any other one. Linking libraries 
of MC-TESTER, ROOT or HepMC is an option convenient also for an interpretation/verification 
of results. 

Basic features which are necessary for our solution have already been tested for Belle and 
BaBar software environments. Changes presented in Appendix [B] are from that perspective 
rather straightforward. 

7 Theoretical basis of the currents 

7.1 Resonance Chiral Theory framework 

The Lorentz structures of currents are universal and a proper minimal set of form factors 
only needs to be chosen for a particular decay mode. One can recall the QCD symmetries 
to gain some insight in the form factors. It is particularly useful that the chiral symmetry 
of massless QCD allows one to develop an effective field theory description, xPT, valid 
for momenta much smaller than the p mass p^ [20] . However, xPT cannot provide pre- 
dictions valid all over the r decay phase space [65] . it constrains nonetheless the form and 
normalization of the form factors in such limiio. 

The computations done within Resonance Chiral Theory {RxT [22l [23] ) are able to 
reproduce the low-energy limit of xPT up to NLO |f2| and demonstrate the right falloff [ 671 
[68j in the high energy region. The current state-of-the-art for the hadronic r decays form 
factors (Fj) is described in Refs. [^ [69]. 

The description provided by RxT complies with the low-energy properties of the under- 
lying theory, at least up to next-to-leading ordeilfj. Nonetheless it is necessary to extend 
xPT to the intermediate energy region, which is probed through hadronic r decays. In 
order to do this, one relies on the large- A^c expansion of QCD ^3\ [TH [75] which predicts 
that in the Nc — ?■ C)0 limit there is an infinite tower of zero-width resonances experiencing 
among them local effective interactions taken at tree level. We model this setting with a 
spectrum that resembles the measured onqf^l- We introduce a NLO effect (whose impact 
is, however, quite sizable in phenomenological applications) in that counting, providing the 



^^One should keep in mind that the effect due to different and nonzero masses of u and d quarks already 
now seems to be necessary to explain some aspects of r data. 

^^The improvement obtained when doing this can be appreciated in Fig. 1 of Ref. |66j . 

^'^Also the leading next-to-next-to- leading (NNLO) order terms [701 [H] are reproduced [5T|[72]. 

•^^In RxT all nine mesons [p, K*, uj, (</>)] have the same mass in the Nc — >■ oo limit without taking 
into account SU{3) breaking. The 1/Nc corrections make the 4> heavier [through the axial anomaly that 
breaks f/(3) down to SU{3)]. SU{5) breaking makes Mp ^ Mr- + M^. 
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resonances with an energy- dependent width computeco within RxT (see Section [3]). 

The appropriate falloff at large energies [671 EH] is imposed to the form factors, vector- 
vector and axial- vector — axial- vector correlators [46l [76], [771 [2HI [79] . This results in a set 
of relations among the coupling constants of the theory, which are obtained working in 
the single resonance approximation (only the lightest multiplet of resonances is included 
per given set of quantum numbers) and in the A*",:; — )■ oo limit. Upon integration of 
the resonances, this procedure allows the saturation of the values of the xPT low-energy 
constants at 0{p^) and O(p^), both in the even- and odd-intrinsic parity sectors. 

The results for all hadronic currents have been calculated within RxT working, with 
the exception of both two-pion and two-kaon modes, in the isospin limit. Therefore the 
corresponding hadronic form factors depend only on the average pion and kaon masses: 

m^ = (mTfO + 2 ■ m^+)/3 , rriK = {jtiko + mK+)/'2. ■ (42) 

For the three-pseudoscalar modes every hadronic form factor consists of 3 parts: a chi- 
ral contribution (direct decay, without production of any intermediate resonance), one- 
resonance and double-resonance mediated processes. The precise form of the form factors 
for three-pion and (two pions - one kaon) modes is presented in Sections 12. ![ 12.21 and 12.31 
For the two-meson modes the corresponding vector form factors are built from the lightest 
resonance contribution in RxT [51]. FSI are resummed by means of an Omnes function 
[51], see Section [231 Both two- and three- meson r decays are sensitive to the exchange of 
excited resonances, whose contribution we have to account for. The exchange of heavier- 
resonances could be computed in the same fashion as for the lightest multiplet, giving rise 
to new unknown couplings^l. This possibility has been postponed for the moment to pre- 
vent the increase in the number of free parameters. We believe that it is sensible to tackle 
this task once more knowledge on the couplings of the RxT is achieved. In the three- 
meson decays, the excited resonances have been included phenomenologically, introducing 
an additional parameter, /3p/, in Eq. ( ITTj) . This has been done in such a way to keep the 
chiral limit result and the QCD-ruled short-distance behaviour. For the two-meson pro- 
cesses, they have been included analogously as the lightest resonance contribution, making 
sure that the appropriate low- and high-energy limits are not spoiled and that unitarity 
and analiticity hold perturbativel}|2j. Following the approximation proposed in Ref. [53], 
three new parameters, related with the p(1450), p(1700) and i^*'(1410) couplings, appear: 
7 = -F(.G"^/F2 [for p(1450)], 5 = -F{:.G'{,/F^ [for p(1700)] and 7x. = -F^G'y/{FFK) 
[for ir*'(1410)] Ej in Eqs. <^^ and d^ZD- The short-distance QCD constraint for the vector 
form factor will require the relation FyGy + FyG'y + FyG'y + ... = F^ to hold. 



^^The width of a spin-one resonance is defined [SS] as tlie imaginary part of the pole generated by 
rcsumniing those diagrams, with an absorptive part in the s-channel, that contribute to the two-point 
function of the corresponding vector current. 

3''See, for instance Refs. [SOlJIl]. 

''^See subsection 17.31 points 1 and 3 and Appendix IeI for an improved treatment. 

^^Due to SU(3) breaking effects 7 7^ "fKn- 
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7.2 The error associated to the 1/Nc expansion 

Any model based on the 1/Nc expansion will naturally raise the question of the error 
associated to that expansion, which can be naively estimated as 1/Nc ~ 30%. It is 
natural to object to the convergence of a series in 1/Nc that ends up being 1/3 in the real 
world. The associated error of the size of the expansion parameter would be much larger 
than the statistical error of experimental measurements or even of our precision target 
of a few percent. Just to give a counterexample, a look at the results of Ref. |5T] will 
quickly suggest that the actual error can be much smaller, which is noteworthy, especially 
taking into account that in the quoted reference, an impressive agreement with data in 
e^e~ — )■ '7r"*"7r~ up to s = 1 GeV^ was obtained in terms of just one parameter, M. 



p- 



Let us recall how this can be possible. First, the large- A''^ expansion of QCD is not 
an expansion in the usual perturbative sense. For instance, in QED, the perturbative 
expansion means, firstly, that the diagrams with less photon couplings to fermionic lines 
dominate. After that, when we compute diagrams both at tree level, and including loops, 
we realize that the expansion parameter is a = e^/(47r)^ and the coefficients of the series 
are small compared to a: every order we go further in the expansion the error reduces 
by ~ 1/a ~ 137. In the large-A"c expansion of QCD, we know first, which diagrams are 
leading-order (planar diagrams with gluon exchanges) and which ones are suppressed. Un- 
fortunately, the expansion has a fundamental subtlety that prevents one from determining 
the expansion parameter after that: there are infinite diagrams at any given order in the 
expansion, so that one cannot perform a calculation at both LO and NLO and compare 
them to know what is the expansion parameter. It is known, however, that diagrams with 
internal quark loops are suppressed as 1/Nc and that non-planar diagrams are suppressed 
as 1/Nq. Moreover, if the number of quark fiavours, rij, is not considered to be smaller 
than Nc the former diagrams can even scale as rif/Nc- These kind of contributions would 
be responsible for the mixing of qq and qqqq states. The fact that this effect is not observed 
in Nature and the success of the quark model classification of mesons in qq multiplets sug- 
gests that the coefficients of these diagrams with internal quark loops are tiny, in such a 
way that the first non-negligible correction would come from the non-planar diagrams, sup- 
pressed as 1/N^ ~ 10%. This reasoning may explain why the large- A"c expansion is a such 
a good approximation for low and intermediate-energy QCD, given the phenomenologi- 
cal successes of its applications in meson effective field theories [82] and the corroborated 
predictions given by the large- A^c limit both for xPT [83l [Hj and for RxT [23] coupling 
constants. All these reasons seem to suggest that, quite generally, some factor comes to 
complement 1/Nc for the value of the expansion parameter to be reduced and the relative 
accuracy to be increased. 

This conclusion is supported by the investigation of some 0{p'^) and O(p^) couphngs of 
xPT, which are done modelling the NLO expansion in 1/Nc of RxT [85], |86l [871 EHl 189] . 
According to the size of the corrections, we judge that 15% can be a reasonable general 
estimate (see, however, our discussion in Appendix O on the possible variations on the 
predictions of the couplings obtained in the Nc — >■ oo limit). Noticeably, the actual 
expansion parameter can be computed for RxT in the study of the vector form factor of 
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the pion at NLO in the 1/Nc expansion [90], yielding 

which, at lowest order, is the ratio of the vector width and mass, ay ~ 0.2, agreeing with 
the previous discussion. 

Moreover, we should emphasize that our approach goes beyond the N^ ^ oo limit. We 
supplement the lowest order in the 1/Nc expansion for the theory in terms of mesons by 
the leading higher-order correction, namely by including the resonance (off-shell) widths for 
the wide states p, K* and Oi. This seems to point to smaller errors than those characteristic 
of the LO contribution in the 1/Nc expansion and may be able to explain, altogether, an 
eventual fine agreement with data. 

7.3 Other sources of error 

Once this major concern on the reliability of the RxT hadronic currents has been discussed, 
let us consider in turn other possible sources of error in the different hadronic r decay modes 
considered. 

1. In the two-meson r decay modes an Omnes type of resummation is employed for the 
FSI. The proposed expressions, Eqs. 1^^ to ( p7|) . respect unitarity and analiticity 
only in a perturbative sense (this will be improved along the lines discussed in Ap- 
pendix [E]). The effect of these violations is, however, pretty small, as one can see 
comparing the results of Refs. [27] and [50] for the r — )■ KttUt decays. That is why 
we consider our current parametrizations for these decay channels a reasonable tem- 
porary approach. For a future update of the program we recall that Eqs. (^^ to ( ITT]) 
should be replaced, using the procedure described in Appendix [E] For the Kir vector 
form factor we are going to uscj the parametrization of Ref . [50] in a future upgrade 
of the program. Analogous works for the '7r~7r*' and K~K^ vector form factors are 
under way. 

2. In the three-meson modes, scalar and pseudoscalar resonance exchange has been ne- 
glected. The spin-one character of the SM couplings of the hadron matrix elements in 
the r decays implies that the form factors for these processes are ruled by vector and 
axial-vector resonances. Their contribution should be minor in r — )■ KKtcUt decays 
(see the related discussion in Section 2 of Ref. [2S]) while that of the scalar resonances 
can be a bit more important in r — ;■ Svri/^ decays. The lightest scalars are, however, 
suppressed in the large- A^c" limit (see, however, Ref. [92]). In these decays, we have 
also neglected systematically three-body FSI for the moment. They may be impor- 
tant in the available phase space and at the required precision [931 [M], thus, should 



We note that another interesting approach, based on Omnes integral equations incorporating both 
chiral constraints at low energies and QCD short-distance relations at high energies was performed in 
Ref. [91]; isospin- violating corrections were also studied. 
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be investigated at a time of comparison to experimental data. The present step of our 
work is devoted predominantly to establish a technical environment. Consequently, 
the effort was concentrated on those theoretical aspects which bring difficulties in 
program design and a careful discussion of potentially less important effects is del- 
egated to the future work. Finally, the computation of Taiiq"^) through the optical 
theorem [21] does not give the corresponding real part of the loop function, which 
we have disregarded. Although this approximation might be supported numerically 
it induces a small violation of analiticity. 

3. In the r — )■ KKtcUt- decays, for the moment, we have not included the contribution 
of excited resonances. We expect that their influence is bigger in this case than in 
the three-pion channels, since the mass of the hadronic system is larger. The impact 
of this error can be as large as the one coming from the 1/Nc expansion. 

7.4 Numerical estimates of the errors in the different decay chan- 
nels and distributions 

What is the numerical precision we should expect in confrontation of our currents with 
the data? As already mentioned, a first crude estimation gave us an error of the order 
of the expansion parameter, l/Nc, thus 30% precision tag for our results, but we already 
accumulated, during the preceding discussion, a number of indications pointing to better 
accuracy of the results obtained within this approach. It is thus of importance ^95j to 
evaluate solely on the basis of theoretical considerations, for which decay channels we 
expect precision to be better or worse and in which regions of the phase space. This is of 
course a must for scientific theory as Resonance Chiral Theory is supposed to be. 

The answer, if scientific theory holds, always comes from the confrontation with the 
data. Results from theory have to be prepared in such a manner that agreement confirms 
and discrepancy invalidates the theory under consideration. Theoretical results have to 
include estimation of their errors from within theory itself. That is the basic principle of 
science methodology and our work has to keep this aspect in mind too. We have to address 
what is input and where such confrontation may take place. 

Until now we see that the precision can be at the level of a few perceno for the two- 
meson modes [6l [271 SHI EOl |53l |96] . Despite the accuracy is at a comparable level for the 
dV /dq^ distributions in the three-pion decays, the error on the distribution in the tt+tt" 
invariant mass is at the level of ~ 20%. We expect to improve it to a few percent level once 
FSI are accounted for. In the r — )■ KKirUr decays the situation is, somehow, reversed. 
The largest error (~ 30%) comes on the dT/dq^ distribution, while the accuracy is much 
better in the K'^tt~ (~ 10%) and K~^K~{^ 5%) distributions. We expect that the errors 
on the KKtc modes are reduced a factor of two when the excited resonances contributions 
will be taken into account. 



^^Once the parameters entering Eq. (j26|) are fitted, we expect similar accuracy for the KK modes as in 
the TTTT and Kt: cases. 
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To get form factors with substantially better agreement with the data, one may need to 
introduce ad-hoc factors, hopefully close to unity and hopefully based on educated guesses. 
This may provide a valuable hint for future theoretical work. 

8 Summary 

In this paper we have documented a set of currents based on Resonance Chiral Theory 
for use in hadronic r decays into either two {tt^tt^, tt^K^, tt^K^ and K~K^) or three 
{ti^ti^ti^, 7i^7[^7c~ , K^TT^K^, K^TT~K^ and K~tt^K^) pseudoscalars. The set covers more 
than 88% of total hadronic r width. Technical tests of the installation into the FORTRAN 
program have been documented. The set can be used as an upgrade easy to install into 
any version of the TAUOLA r decay library. In this way currents are ready for confronting 
the r decay data (unfolded or not). Precision fits can be performed and arrangements for 
use of model-dependent weights are ready for that purpose. 

On the technical side, the r decay algorithms themselves have been checked down to 
0.05% precision level. To this end, a detailed comparison between analytic and Monte Carlo 
results has been provided. Statistical samples, 2 (3) orders of magnitude larger than at time 
of reference [H] ([H]) were used. This technical precision of 0.05% is substantially better 
than physics precision of our currents which we estimate at the 5-30% level, depending on 
the channel. Software environment for further phenomenological work is prepared. 

In the present work, we have concentrated on decay modes contributing to 88% of the r 
hadronic decay width. In particular, the decay modes of r into 47r's have not been updated 
in our work. The currents for these decays, constructed on the basis of low-energy e~^e~ -^ 
hadrons data have been available since Ref. [33]. For the time being this can be used as 
an alternative to TAUOLA cleo default of Ref. [14] in discussion of the systematic error for 
background modelling. In these decay modes theoretical foundations are less profound, 
but we hope that the presented solution will not jeopardize analysis of other channels due 
to lesser control of cross contaminating channels. Together with 47r's currents described in 
Ref. p3] our system covers now 97% of the total hadronic r width. 

It is straightforward to extend our work to more elaborated currents including, e.g. 
scalar form factors. Because of different nature of benchmark distributions, technical 
aspects and different physics assumptions we leave this task to the forthcoming work ^22] . 

As a consequence of comparisons with the data, some of the theoretical assumptions 
may need to be reconsidered too. We start from a theoretical approach common for all 
channels. Empiric form factors may need to be added later though. In this way not only 
agreement with the data will be established, but eventual inefficiencies of our starting 
approach will be numericaly evaluated. That is why, it is important to ensure that model- 
dependent weights can be calculated after (and independently of) the detector effects. 
There are several assumptions which lie behind our calculations and in case of discrepancy 
with the data they may need to be revisited. Let us list them in descending order of their 
theoretical foundation: 

1. We assume that Lorentz invariance will not need to be reconsidered at any step of 
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our project. 

2. With respect to the separation of the matrix element into leptonic and hadronic 
parts, we assume that electroweak corrections will not affect such separation beyond 
precision level of several permille at most, the question of an overall normalization 
factor for all hadronic channel, see e.g. Refs. [971 [981 [99] , is of no practical importance 
for Monte Carlo. 

3. We assume that isospin symmetry should be a good guiding principle. We suppose 
that it should hold more accurately for distributions than for amplitude phases, but 
we do not expect large effects. 

4. Some of the effective couplings can be predicted by considering the asymptotic be- 
haviour of Green functions and form factors both in the effective theory (RxT) and in 
the operator product expansion of QCD [23j. However, these predictions are affected 
by different sources of errors, most importantly the model dependence on the real- 
ization of the large-iVc limit of QCD [TOQl [Ml [M [I03l ^QM EQB [1061 [M] (mainly 
the choice of the resonance spectra, but not only). Special care should be taken when 
relating different channels, especially if the statistics in both of them is very different 
and in one decay channel only a subset of resonance parameters is used. One should 
not forget that the (formal or not) integration of heavy degrees of freedom out of 
the action affects the values of the parameters in the remaining lower-energy theory. 
Results for the individual modes should be analyzed consequently. 

5. Our effective couplings and interactions are based on the low-energy effective field 
theory of QCD (xPT), whose results are reproduced at NLO in the corresponding 
limit by RxT. Although the latter being formally sound, there is model dependence 
in any realization of the large- A'',^ limit of QCD for mesons and, moreover, we are 
introducing the contribution of excited resonances only at a phenomenological level, 
see Eq. ( 1TT1) and the related discussion in Section [71 a feature that can be improved 
in the future. We may need to explore the limits of such approach and take feed back 
from experiments. 

6. Different solutions have been advocated for taking unitarity properties, via the prop- 
agator widths, into account. In particular, there is no consensus that the exponenti- 
ation of the real part of the resonance width is always the best solution. This point 
must be investigated further, especially in light of precision fits to high-statistics 
data. 

At the same time, the seeming violation of these principles may be a consequence of 
some experimental problems. That is why such discussion requires simultaneous partici- 
pation of theorists and experimental physicists and the proper software environment. 

Such fine-tunings took place for TAUOLA version of Ref . [TT] . The improved agreement 
with the data of CLEO and ALEPH has been achieved thanks to the effort of these col- 
laborations resulting with TAUOLA cleo and TAUOLA aleph initializations. The code is 
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available from Ref . [H] . We should expect similar work with our present parametrizations. 
We hope that, in this way, theoretical and experimental constraints will be appropriately 
matched leading to a representation of experimental and theoretical achievements in com- 
mon language of value for future research. Once such confrontation with the Belle or 
BaBar data is completed, one will profit from the technical precision of the simulation 
established here and thanks to prepared flexibility for FORTRAN [H] and C++ [62] users of 
high-energy experiments such as at LHC or for LC will benefit as well. If analytic form of 
the form factors will be attained, it may help start future work for such research as lattice 
calculations. 

We consider this work as a step towards a theoretically rigorous description of hadronic 
tau decay data but, at this moment, we can still not be sure that the currents for all 
channels will be able, after fits, to describe all data well. In this case, detailed numerical 
information on the offending distributions will be provided. Further theoretical work is 
stimulated. Some intermediate results are presented in |1U8[ I1U9| IllUj . 
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A Useful functions and notations 

To ininiinize repetition and to reduce the size of formulas in Section [2] the lengthy ones 
were moved to this appendix. 

In the description of the three-hadron currents the following functions were used: 

A''{q^,x,y,ml,ml,ml) = 3x + ml - ml + (l - ^;^j [2q^ - 2x - y + ml - ml] , 

B^{x,y,m\,ml) = 2 [ml - mj) + (l - ^;^j [y - x + ml - ml] , 
A''''{q^,x,y,ml,ml,ml) = {X' + X") {-3 x + ml - ml) 

(Z 771 \ 

C^{q'^,x,ml,ml,ml) = {ci - C2 + c^) q"^ - {ci - C2 - c^ + 2cq) x 



(ci + C2 + 8c3 - Cs) ml + 8c4 {ml - ml) , 



C™(g^x,m^) = 4 (g^ + x) + (c/i + 8 c/a - 4) m^ 



D'^iq^x.y) = {gi + 2g2-gz){x + y)-2g2{q^ + m]^) 

-{gi - 93) (3 mj^ + ml) + 2 g^ {m\ + m^) + 2 g^ m\ , 

E"(x,y) = {g, + 2g2-g^){x-y). (44) 

They follow conventions of [25]. Function H{x,y) is defined in subsection 12. ![ formula ([7]). 
In the description of the two pseudoscalar form factors, following [121 120] , the function 
Apq{s) was used to describe loops involving pions, kaons and 1] mesons 

^^^l^j__192^VW£)^W£W, (45) 



where 



1 , .„..,. aIq 1 



Mpq{s) = -^{s - 2^pq)Jpq{s) + -:f^JpQ{s) - ^kpg + -^ , (46) 



F2 



j:pQ = mp + mQ, ApQ = mp-mQ, ^pq = {^p - fig) (47) 



and 



ApQ 

Lpq{s) = ^f^JpQis) . (4J 
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Finally, 



o2 



The /ip = 32^2^a In ( —& ) (at present we take /i = Mp for all 2 pseudo-scalars mode: 



1 



ApQ J:pq{s)\^ ml u {s + uf-AlQ 



, In — ^ In ■ ,„ ^ „ 

ApQ y m|, s [s- uy - Af,Q 



327r2 
Jpq{s) = Jpq{s)-sJ'{0), (49) 

where z/^ = A(s,mp,mQ). Care has to be taken to keep the imaginary part of z/^ in the 
phase-space regions where u"^ < 0, and where Im v must not be set to zero, see the function 
JPQ1_FUNCT in the file funct_rpt.f . In the vector form factor for two pions, this effect 
shows up in the KK contribution to the loop function from below the threshold of KK 
production. 

For the p' and p" mesons the loop function A^^i^s) taken from [51] reads: 



11^ ) s 3 \a. 



AM = In ^ + 8^ - - + a;J In ^-— . (50) 



In the last formulas the SU{2) limit is taken and m^± = m^o = 0.13804 GeV. 

All functions in this appendix are coded in file new-currents/RChL-currents/f unct_rpt . f 

B Installation 

Our project tar-ball, even though resulting from a rather large effort, is not designed for 
independent installation. This would be of course straightforward and we will return to 
that solution in the future, once the currents are optimized to improve agreement with 
the data. The parametrization will become integrated part of the TAUOLA distribution; for 
fortran or for C++ use, like in references [T41IB2] . At present we concentrate on a solution 
which is most convenient for the experimental user e.g. from Belle or BaBar collaboration 
aiming at combining the code with the version of TAUOLA which is already being used as 
part of the simulation set-up. We aim at preparing an add-ujo for already existing set-up. 
The tar-ball can be downloaded from the Web page [28J of our project. 

Once tar-ball is unpacked inside TAUGLA-FORTRAN/tauola subdirectory (of user envi- 
ronment), the directory tauola/new-currents will be created, all necessary fortran files 
will be found there. For convenience, later on, we will use the following aliases: 

• ${RCHLCURRENTS} instead of tauola/new-currents/RChL-currents . 

• ${GTHERCURRENTS} instead of tauola/new-currents/other-currents . 

• ${ INSTALLATION} instead of tauola/new-currents/Installation . 



■^^ Alternatively, for n we raay take Mp for PQ = tttt, KK and Mk' for PQ = Kit [2^. 
^■^As the project is developed under svn, the tar-ball is accompanied with svn label and it should be 
kept for reference. 
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In tauola/new-currents further sub-directories for more advanced use or for docu- 
mentation will be found: 

• ${RCHLCURRENTS}/tabler/al - programs for pretabulations in particular of g^-dependent 
ai widtlQ. 

• ${RCHLCURRENTS}/cross-check - code for technical and numerical tests. 

• new-currents/paper - present paper. 

• ${ installation} - instructions for modifications to be introduced in FORTRAN files 
and makefile residing in directory tauola . 

• new-currents/Installation-Reweight - instruction and example of using reweight- 
ing algorithm. 

None of the directories listed above contains code which is to be loaded together 
with TAUOLA library. Code loaded with the library is located only in the main folder 
of ${RCHLCURRENTS} and tauola/new-currents/other-currents . Programs in 
${RCHLCURRENTS}/tabler/al can update the fortran code located in file 
${RCHLCURRENTS}/initAlTab.f . 

Once installation is completed, to invoke the calculation of our new currents the CALL 
INIRChL(l) has to be invokecoby user main program prior to call on TAUOLA initialization. 
If instead CALL INIRChL(O) is executed prioiQ initialization, old currents - as in Ref. [H] 
- will be used in generation. 

The CALL INIRChL(l) may activate also new currentO e.g. for r]nn or Att decay 
channels. At present only wrappers of currents of Ref. [H] and |33] are prepared in the 
directory ${OTHERCURRENTS} for convenience of users and our future work. These currents 
lead to substantially different distributions, that is why, one may require adjustment of 
phase space presampler used to optimize speed of generation. Anyway, as a default, they 
are turned off. For INIRChL(l) the same TAUOLA cleo currents as for INIRChL(l) are 
used. To turn other options, ISWITCH located in the file ${OTHERCURRENTS}/f f ourpi.f 
has to be changed from its default value to 1, 2, 3 or 4. 

B.l Changes for host TAUOLA version 

In order to use new currents, changes have to be made to the host TAUOLA installation. 
Let us document here in great detail changes to be introduced in TAUOLA cleo version. If 



^'^In the directory ${RCHLCURRENTS}/tabler the place to calculate other pretabulated functions (as 
possibly the scalar form factor for Kn channel) is reserved. 

''^For a C++ user, examples of use of inirchl_(l) are given in 

new- current s/ Inst allat ion-Rewe ight / directory. 

■^^ This can be done also after initialization as no initialization of tables is needed. Then one can revert 
the change again with CALL INIRChL(O) . 

"^^ Although the t — > r^^'V^Tr^i^T decays have been worked out within Resonance Chiral Theory |111] . 
the corresponding expressions for the currents have not been incorporated yet to the program. 
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some modifications were introduced and user's host TAUOLA installation differs from TAUOLA 
cleo of Ref. [H], then modifications prepared in ${INSTALLATI0N} directory can not be 
used directly and some adaptation may be necessary. In either case we advice to check if at 
least some of the numerical results from Ref. |28] are correctly reproduced after installation. 
Let us list now changes which have to be introduced to files residing in TAUOLA/tauola 
directory of the user installation . 

• TAUOLA/tauola/makefile 

The list of LIB_OBJECTS must be extended and additional objects added: 
${RCHLCURRENTS}/f 3pi_rcht . o , ${RCHLCURRENTS}/f kkpi . o , 
${RCHLCURRENTS}/f kkOpiO . o , ${RCHLCURRENTS}/wid_al_f it . o , 
${RCHLCURRENTS}/f rho_pi . o , ${RCHLCURRENTS}/f unct_rpt . o , 
${RCHLCURRENTS}/value_parajneter . o , ${RCHLCURRENTS}/initAlTab . o , 
${RCHLCURRENTS}/f kpipl . o , ${RCHLCURRENTS}/f kOk . o , 
${OTHERCURRENTS}/f etapipi . o , ${OTHERCURRENTS}/f f ourpi . o , 
${OTHERCURRENTS}/binp . o , ${OTHERCURRENTS}/curr_karls . o , 
${OTHERCURRENTS}/curr_karls_extracted . o; 

if there are no additional dependencies the ${lNSTALLATION}/makef ile-tauola file 
can be simply copied into tauola/makef ile . 

• TAUOLA/tauola/tauola.f 

If the file in user's version coincides with the one of TAUOLA cleo distribution, the 
${lNSTALLATION}/tauola.f-new file can be simply copied into tauola/tauola.f . 
To verify this, the diff file ${lNSTALLATION}/tauola.f-oldDIFFupdated may be 
inspected. 

• TAUOLA/tauola/f ormf . f 

If the file in user's version coincides with the one of TAUOLA cleo distribution, the 
${lNSTALLATION}/f ormf .f -new file can be simply copied into tauola/f ormf .f . 
To verify this, the diff file ${lNSTALLATION}/f ormf .f-oldDIFFupdated may be in- 
spected. 

Once changes are introduced the new currents will be activated and the old ones will be 
overruled once call to routine INIRChL(l) is invoked. Otherwise, or if CALL INIRChL(O) is 
invoked (at any time), old currents will be then switched back on. The routine INIRChL(l) 
has to be invoked by the user program at the initialization step. For the C++ user, a 
definition of extern "C" void inirchl_(int i) ; has to be included and execution of 
inirchl_(&i) ; performed. 

An example has been provided in ${lNSTALLATION}/demo-standalone . It is based 
on default TAUOLA cleo example with the only modification being the call to INIRChL(l) 
before default TAUOLA initialization. 
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B.2 Calculating numerical tables used by form factors 

The directory ${RCHLCURRENTS}/tabler/al contains the program dalwid_tot_rhol_gauss .f ; 
it creates a table of ra^(g^) according to Eq. (!35ll . The system does not use any information 
from TAUOLA initiahzation except the pion and kaon masses. 

We have also prepared a place to add tables for other functions, in the near future it 
will be done for the scalar form factor of the Ktt mode. 

B.2.1 Executing the code 

The program dalwid_tot_rhol_gauss.f produces the q^ distribution of the ai off-shell 
width. 

To compile, type make in ${RCHLCURRENTS}/tabler/al directory. To run, type make 
run. Each hne of the produced output includes the value of g^[GeV^], and the value of 
dT/dq^[GeY~\ This table is written into the file initAlTab.f which is the FORTRAN 
code ready to use. One can shift it to ${RCHLCURRENTS} directory by make move command. 
Text format table is written into file widal_qq_tot_2e5 . out . 

B.2. 2 Setup 

Input parameters and common blocks are located in ${RCHLCURRENTS}/parameter . inc . 
Other parameters are defined in ${RCHLCURRENTS}/value_parameter .f . These param- 
eters may be changed by the user. If the parameters affect the g^-dependent ai width 
(or other pretabulated functions), the tables need to be generated anew with the help 
of programs residing in the directory ${RCHLCURRENTS}/tabler . A list of the param- 
eters that affect generated tables (and thus require tables to be generated again) is in 
${RCHLCURRENTS}/value_parameter .f . Some of the variables used in functions from 
${RCHLCURRENTS}/funct_rpt.f are declared in ${RCHLCURRENTS}/funct_declar. inc . 

B.3 Tests 

Directory ${RCHLCURRENTS}/cross-check contains three subdirectories: 

1. check_analyticity_and_numer_integr, it includes: 

• test of numerical stability in calculations of r^^ (g^) and for the whole r hadronic 
decays as described in Section [2l For that purpose it is checked if continuity of 
results as a function of the invariant mass holds. 



ri 



the result for the integrated width of the r — )■ 27^/,-, r — )■ KttUt-, t -^ K^K v. 
T — )■ Stti/t-, r — )■ Kn'Rur and r -^ K'TT^K^Ur- These results can be con- 
fronted with the result of Monte Carlo simulation collected in sub-directory 
tauola_result_modes . 
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2. results_numer_iiitegr_3pion presents the results for the width of r — )■ Svrz/^ as a 
function of the 3 pion invariant mass. It is calculated by numerical integration of the 
analytical formula for different choices of hadronic form factors as it is described in 
Section HI 

3. tauola_result_modes contains Monte Carlo results for both differential and total 
width for the processes r — )• 27ri/^, r — )■ KnuT-, r — )■ K'K^Ut-, t — )■ StiUt-, t — )■ 
Kn-Kur and r -^ R-ti^K^u^. 

B.3.1 Numerical stability tests 

The directory ${RCHLCURRENTS}/cross-check/check_analyticity_and_numer_integr con- 
tains six subdirectories with tests of numerical stability for hadronic r decay modes and a 
subdirectory with the test for the ai width. Each decay channel is located in a separate 
directory. Details regarding each of these tests are described in README files of the directory 
and every subdirectory as well. That is why only the basic information is provided in our 
paper. We have checked using interpolation from neighbouring values that the value of 
dV /dq^ is continuous and is not contaminated by numerical instability of multidimensional 
Gaussian integration. Also we present the analytical results for the partial width of every 
channel to be compared with the Monte Carlo ones. 
Content of the directory: 

• check_analyt_3piwidth: test of numerical stability of the distribution dT{T — )■ 
z/T-7r7r7r)/c/g^. Results are presented for separate modes: dVij — )■ z/T-7r°7r°7r~)/dg^ and 
(ir(r —7- z/^7r^7r^7r+)/(ig^. Also the value of the partial widths for the channels is 
provided for the comparison with the TAUOLA results. 

• check_analyt_kkpi - test of numerical stability for dT{T — t- v^KK'K)/dq^. Re- 
sults are presented for separate modes: dVij — t- v^K~ ti^ K^) / dq^ and dVij — ?► 
Vt-K^t[~ K^) / dq^ . The value of the partial widths for both channels are provided. 

• check_analyt_kkOpiO - tests of numerical stability for r — )■ Ut-K^tt^K^: both the 
spectrum dr{T — )■ Vt-K~ tt^ K^) / dq^ and the partial width are provided. 

• check_analyt_2pi - tests of numerical stability for r — )■ z/t-tt^tt'^: both the spectrum 
(ir(r —7- UrTT' TT^) / dq"^ and the partial width are provided. 

• check_analyt_kpi - tests of numerical stability for r — )► UrKn: both the spectrum 
for the total width (ir(r — )■ VrKu) / dq^ and the partial width for channels it~K^ and 
TT^K~ are provided. The partial widths for the individual decays are checked to be 
2/3 and 1/3 of the total Ktx width, mass effects are negligible in this case. 

• check_analyt_kOk - tests of numerical stability for r — )■ UrK'K^ both the differential 
distribution (ir(r — )■ VrK~ K^) / dq^ and the partial width are provided. 

• check_analyt_altable - tests of numerical stability of Ta^{q'^) produced by program 
described in Appendix IB.2I 
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B.3.2 Analytic integration test 

The results of the analytical integration test in the three-pion case are presented in the 
directory ${RCHLCURRENTS}/cross-check/results_nunier_integr_3pion . They are pro- 
duced by the program totwid3pi_qq_table.f in the directory 

${RCHLCURRENTS}/cross-check/check_analyticity_and_numer_integr/check_analyt_3pi . 
The program can be compiled by command make and run with make totwid3pi_run > 
output . txt . 

The setup file input_flf2f4.dat, in 
${RCHLCURRENTS}/cross-check/check_analyticity_and_numer_integr/check_analyt_3pi 
contains: 

• eps - defines (relative) precision of the Gaussian integration. 

• kf 1 - fiag for form factor Fi. For kf 1= 0, 1 or 2 Fi will be set respectively to 0, 1 
or to its functional form. 

• kf 2 - fiag for form factor F2. For kf 2= 0, 1 or 2 F2 will be set respectively to 0, 1 
or to its functional form. 

• kf4 - fiag for form factor F4. For kf 4= 0, 1 or 2 F4 will be set respectively to 0, 1 
or to its functional form. 

• Chan - fiag to choose the 3 pion mode. chan= 1 for tt^tt^tt" and chan= 2 for 7r~7r~7r"'". 

If the functional form of the form factors is used, it will be taken from the file 
${RCHLCURRENTS}/cross-check/check_analyticity/check_analyt_3piwidth/funct_3pi .f 
If kf 1 or kf 2 is set to 2, pretabulated file ${RCHLCURRENTS}/initAlTab.f will be used for 
Fai in the propagator of the ai-mesoro. 

Output file contains four columns: 

• qmin (in [GeV^]) - lower boundary for the integration over 3-pion invariant mass. 

• qmax (in [GeV^]) - upper boundary for the integration over 3-pion invariant mass. 

• eps - estimate of the integration precision in the result. 

• total width (in [GeV]). 

Results for the different configurations of the form factor are presented in 
tauola/RChL-currents/cross-check/results_numer_integr_3pion . 



*Note that the tabulated file is generated by the program described in Appendix IB. 2 1 
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B.4 TAUOLA weight recalculation mode 

Let us present now the installation necessary for the method of weighted events, which 
was envisaged in Section 16.11 An example of such installation code is included in our 
distribution tar-ball in directory new-currents/Installation-Reweight . 

Before reweighting method can be used, TAUOLA needs to be adapted to new currents as 
explained in the Appendix lB.il Afterwards, our example program tau-reweight-test-ASCII . c, 
residing in the directory new-currents/Installation-Reweight, can be run with the help 
of the simple make commancO- For more details regarding the reweighting examples, refer 
to README located in new-currents/Installation-Reweight . 

The following subsection describes reweighting algorithm as well as initialization used 
in the example. Note, contrary to the rest of the project, reweighting algorithm, including 
examples of its usage, is written in C++. 

B.4.1 Weight recalculation algorithm 

In order to use recalculation mode, several steps have to be performed from user program: 

1. Before TAUOLA initialization, RChL currents have to be switched on. This can be 
done with the help of the wrapper for FORTRAN function INIRChL(IVER), by calling 
inirchl_(&i) ; with i = 1; . Two versions of currents will be used, but initialization 
must be done for IVER=1, for initialization of RChL-specific variables and tables. 

2. Initialization of TAUOLA must be called. We are using initialization taken from the de- 
fault TAUOLA example, stored in wrapper function f _interf ace_tauolaInitialize. 

3. For each event, the information about r and its decay products must be filled and 
stored in instances of SimpleParticle clasqfj. 

4. Once the kinematical configuration for r decay is read from the datafile (or fifo pipe), 
function: double calculateWeight (SimpleParticle &tau, 
vector<SimpleParticle> &tau_daughters) can be used to retrieve the weight. 



The algorithm of the function calculateWeight is sketched in the following: 

1. Particles are prepared and boosted to the appropriate frame. 

2. TAUOLA decay channel is identified. 



^^Other example, tau-reweight-test-HepMC.c, requires installation of HepMC and optional installation 
of MC-TESTER, and their paths provided in the Makefile . 

^^Class SimpleParticle is used only to contain four-vector and flavour of the particle. 
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3. TAUOLA cleo currents are switched on with inirchl_(&i); i = 0. 

4. Call to appropriate internal TAUOLA FORTRAN subroutine, returning weight WTl. 

5. RChL currents are switched on with inirchl_(&i) ; i = 1. 

6. Call to the same routine as in step 4 is performed, returning weight WT2. 

7. Ratio of weights calculated at steps 4 and 6 gives required model replacing weight. 

8. WT = WT2/WT1 is returned to the main user program. 

It is rather straightforward to extend this method to the case when more than one new 
version of physics initialization is to be used. Note that the examples are set up so that the 
weight is calculated both for r~ and r"*" and stored in variables WT_M and WT_P respectively. 
In cases where only a single r is present in the event, the weight corresponding to the 
second r equals 1.0. 

Alternatively, in cases when this approach cannot be used or is inconvenient, variants 
of the method, based on fifo pipes can be useful as well. Prototypes for such solutions can 
be obtained from Ref. |112j . 

Hadronic currents for r+ and r^ differ due to CP parity. The resulting effects are taken 
into account in the reweight algorithm. 

B.5 TAUOLA++ installation 

Thanks to the modular construction of TAUOLA C++ Interface [62], new currents can 
be used in C++ projects in a straightforward way. It is enough to replace the previous 
TAUOLA-FORTRAN installation with the new one, adjusting Makefile with a list of the newly 
added object files. 

For step-by-step instructions, we refer to ${ INSTALLATION} /README-TAUOLA++ . Our 
package has already been tested to work with TAUOLA C++ Interface vl.0.5, but the 
installation procedure is similar for all previous versions and should remain unchanged for 
future versions as well. 

C Input parameters 

The results collected in this paper represent a technical test of program installation as 
well. Figures should be reproduced if the input parameters, collected in Tables [3l HI 
and|5]and defined in routine tauola/new-currents/RChL-currents/value_paraineter .f 
remain unmodified. In some cases the actual numerical value of parameters depends on 
chosen decay channel. For the KKtt modes we use Mk* = {Mk*± + Mk*o)/2. For the 
Kn modes the value of parameters depends on parametrization. It is distinct for the 
one of Ref. [2Z] and of Ref. [5U]. In the second case the mass parameters are noticeably 
different from the pole values. The results for the latter are consistent in both approaches. 
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The choice between the two parametrizations for Kti modes (channels 2, 3 in Table [6]) 
is controlled by FFKPIVEC again to be set in value_parameter . f . The FFKPIVEC = 
activates parametrization from Ref. [50] and FFKPIVEC = 1 the ones of J27j. Numerical 
values of all parameters, not only the masses, are different for the two cases. Variables are 
named with big greek letters for FFKPIVEC = 1 and with the small ones for FFKPIVEC = 
0. 

There are two other flags FFVEC and FFKKVEC in value_paraiiieter .f . The first one 
fixes run with/without FSI effects (FFVEC = 1 for run with FSl effects) and the last one 
chooses the parametrization for two-kaon form factor with/without the excited p meson 
states (FFKKVEC = 1 for the parametrization with p' and p" ). By default FFVEC = 1, 
FFKPIVEC = 1 and FFKKVEC = 0. 

On technical side, the choice of the internal flag KAK is made at the start of each r decay 
generation. It depends on the decay channel labelled by imode (generated by TAUOLA) and 
the flagcEi] FFVEC, FFKKVEC and FFKPIVEC. The variable KAK is then passed into routine 
value_parameter .f and the appropriate choice for the parameters is made. The KAK 
parameter coincides with imode for all channels except the Ktt modes. For the Ktt decay 
modes KAK = 70 if FFKPIVEC = and KAK = 71 if FFKPIVEC = 1. 

For FFKPIVEC = (that is for KAK = 70) parameters marked in Tables with f are used, 
otherwise defaults of Tables H] and [5] are left unmodified. 

For KAK=4 (i.e., for r — >■ 7r~7i^i>r) masses and widths of p,p' and p" result from the 
adjustment to the experimental data and do not coincide with PDG defaults. For other 
channels we simply take the PDG values [3Tj for the p{p') parameters. The PDG values 
are also taken for the narrow width resonances u and 0, numerical values are collected in 
Table H 

The PDG value is taken for the ai masa^^l. The parameters of the Resonance Chiral 
Theory are given in the Table H] as welO. 

The parameters By and Fk can be varied by the user starting from the code version of 
the year 2012. We follow Ref. |2^ and the case of ideal mixing {6v = 35.26°). In this case 
the (j) contribution to r — )■ KKt^Vt- vanishes [25] . However, the absence of intermediate (f) 
exchange contradicts the results of the BaBar Collaboration [113] for the isospin related 
decay e'^e'^ — t- K^ir^K^ and for r decays themselves (see^"'^. The parameter F^ is not 
used in our default formulas. It enters the non-default parametrization for Kit vector 
form factor, i.e., for FFKPIVEC = 0, KAK = 70. We follow Ref. [5U] in the choice 
Fk = 1.198 ■ F. Fk ^ F is related to SU{3) breaking and higher-order chiral corrections. 
The parameter Hto = —1.24004 ■ 10~^ does not appear in the text either. It corresponds 
to the value of the Kn loop function at zero-momentum transfer, Hktt{0), in Eq-(ll) of 



^-'^At the moment KAK depends on FFKPIVEC only. However, in future it will depend also on FFVEC and 
FFKVEC. 

^^For discussion on the difference between the mass used in the resonance Lagrangian and the physical 
one and a possibility to substitute the first with the latter, see [25| and footnote'^ in this paper. 

^•^We point out that the values for the parameters Fp/, Mp», Fp//, 7, 6, 0i and (j)2 lie outside the educated 
guess for its range of variation given in Table ID This is irrelevant for the technical check we are proposing 
in this Section but matters for the actual use of the program. 
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Ref. EOl. 



Parameter Var. name Default Used in channel 



rrir 


MTAU 


1.777 


all,* 


m^r 


MNUTA 


0.001 


all,* 


COs6'cabibbo 


set in TAUOLA init. 


0.975 


all 


Gf 


set in TAUOLA init. 


1.166375-10-5 


all 


m^± 


mpic 


0.13957018 


all,* 


m^o 


mpiz 


0.1349766 


all,* 


TTlri 


met a 


0.547 


2,3,5-9,* 


mK± 


mkc 


0.493677 


all,* 


rriKO 


mkz 


0.497648 


all,* 


M^ 


mom 


0.78194 


7,8 


r. 


gom 


0.00843 


7,8 


M^ 


mphi 


1.019 


7,8 


r^ 


gphi 


0.0042 


7,8 



Table 3: Initialization parameters defined in TAUOLA main code or in file 

new-currents/RChL-currents/value_parameter .f : constants and defaults. In this 
table our defaults used for plots or parameters not requiring to be changed 
in fits are collected. Channels identification numbers are defined in Table 

ini Energy units are powers of GeV. Variables requiring rerun of pretabulation 
new-currents/RChL-currents/tabler/al/dalwid_tot_rhol_gauss . f are marked with *. 

Let us stress that in practice the parameters may need to be varied. The defaults and 
the expected variation ranges are given in Tables S] and O 



C.l Range of variation of the non-resonance input parameters 

As long as the PDG values do not change, the values listed in Table [3] should remain 

unchangeco- 

In order to account for the uncertainty given by higher-order chiral corrections, we 
suggest to vary Fk as indicated in Table HI 



C.2 Range of variation of the resonance input parameters 

The resonance parameters are of different nature in this respect. Apart from the safe 
identification My = Mp J116j . there is more uncertainty and model dependence on them. 
For the program user this is translated in a relative freedom to change the values of Mp, 



^''^The PDG limit on niy^ (18.2 MeV) is not used in the program. If one wants to play with this limit 
[11411115] . our test results shall change in a rather negligible way. We use rrii,^ — 0.01 GeV. 
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Mai, ^P'^ Tp/, Mpn^ Vpii, 7, 5, 01, 02, M^^.i, M^*o, Mi^., M/^./, mi^., m^*/, 7^^*, 7^./, 
r^,', Fy, Gy = F"^ / Fy (although some deviations to this relation -below 20%- may be 
expected due to the effect of excited resonances), F4, /3p and 7_ftv. The changes of these 
parameters can be guided by the educated guesses on their rangCj, displayed in Table HI 
The warning is that the Fy and Fa cannot be changed independently since they should 
satisfy, to a reasonable accuracjo, the first and second Weinberg sum rules taken in the 
single resonance approximation |119] : 

F^-Fl = F\ F^Ml = FlMl . (51) 

Violations of these relations can be due to the modelization of the resonance spectrum in 
the \aige-Nc limit but should remain below 20%. 

There is more uncertainty on the couplings belonging to the odd-intrinsic parity sector, 
namely the Cj, di and gi values and variation ranges are given in Table [51 Some remarks 
on the relations are in place: 

• Ci — C2 + C5 7^ would violate maximally the short-distance QCD-ruled behaviour for 
the vector- vector correlator ny(g^) [25], this condition must not be changed. 



2(74 + g^ comes from r(ci; — )■ tt+tt^tt^). Both the direct production mechanism 
and the one-resonance exchange [16] were taken into account consistently and the 
error is under control: 2(^4 + g^ = —0.60 ± 0.02. If the PDG value for this decay 
width does not change, the value of this combination of couplings should be changed 
within the quoted error only. 

The high-energy large- A^c predictions for the set {5^2, 5'i — 5'3, ci — C2 — C5 -|- 2c6, d^} 
come at the same order in the expansion of ny(g^) in powers of 1/g^ [25]. Therefore, 
changes on the values of these parameters shall be expected from subleading correc- 
tions in 1/Nc and, moreover, they will be highly correlated. Variations of 1/3 with 
respect to the values of any of them, see Table may occur. 

The predictions for di + 8^2 — c?3 and Ci + C2 — 8C3 — C5 were not obtained in hadronic 
r decays but in the study of the (VVP) octet [16] and singlet |120j Green function. 
Therefore, an educated conservative guess yields a deviation up to some 50% of the 
former. For the latter, a non-zero value may arise provided it has a minor effect in 
the observables, since its contribution vanishes in the chiral limit. 



^^Keep in mind, however, the warning concerning the relation Gy — F'^/Fy and the one affecting 
Eq. ([5T|) for the range for Gy and Fa, respectively. 

^^Checking the first of Eqs. ([5T|) is straightforward; for the second one, it should be observed that 
in the different relations among couplings which can be obtained from short-distance QCD constraints 
[67l [68l 11171 1118] and involving Ma, the identification Ma ~ Ma^ is not appropriate |116) . There is 
some tension on the value of Ma: 998(49) MeV in Ref. [81] versus 920(20) MeV in Ref. [85]. The 
range [900, 1050] MeV should accommodate reasonable variations of this parameter in order to estimate 
the possible violations of the second Weinberg sum rule. The interval given for Ma^ is only marginally 
consistent with the PDG value [3T]. However, this is not an issue, since it depends strongly on the precise 
definition of the resonance mass used to extract it; the PDG one and the one in Ref. [24] are different. 
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• Our current understanding seems to point to different values of C4 and g^ tlian those 
shown in Table [51 We would suggest that they are allowed to vary freely in the fits 
(keeping control on the branching ratio if it is not included as a data point in the 
fit). 

In addition to what is explained above, the values of the couplings can be affected by the 
introduction of the second multiplet of resonances in r — )• KKni/r decays. 

Finally there are parameters in value_parameter .f which are already prepared for the 
future update, but not yet used in the program: gro, mf2, gf2, mfO, gfO, msg, gsg. 

D Benchmark results 

In Table [6] we collect results as coming from our default setting. It documents directly the 
distribution tar-ball. Results have to be checked once tar-ball is unpacked and installed 
in the particular environment. It will check if sufficiently large samples are installed, 
correctness of coupling constant setting, etc. This table will be updated in the future, once 
currents are modified, and included in the paper; version of the distribution tar-ball. 

A wealth of data is available from the project Web page [28]. Many of the results 
collected there were obtained using MC-TESTER [SZ]. The distributions can be thus easily 
used for benchmarking any other Monte Carlo program independently if it is used for 
simulations coded in C-I--I- or FORTRAN, provided that event record such as HEPEVT or 
HepMC [122] is used. 

E Final state interactions 

Final state interactions have been taken into account for the two-meson r decays but 
not for the three-meson channels. It is of great importance to include them in the two- 
meson decays, because otherwise the right normalization at the peak is clearly lost and 
the curve would go systematically below the data in the resonance regional. This effect is 
certainly much smaller for the three-meson r decays, where nonetheless it can show up in 
the two-particle invariant mass distributions in particular regions of phase space, where our 
calculation is particularly far from experimental data, see Fig. |H1 For the two-meson modes 
(r — >■ PQur), FSI are taken into account through the Omnes resummation [123] provided by 
the exponents in Eqs. (!24|) and ([27]) . In this approximation [51], the imaginary part of the 
loop function is kept in the denominator, providing the width of the exchanged resonance. 
The real part is resummed in the exponential -while the exponential is a common factor in 
Eq. f l27p . it is different in every term in Eq. f l24p . If we had good theoretical knowledge of 
the energy-dependent meson widths (also for excited resonances) the latter procedure would 
be preferable. For the time being, both approaches are equivalent within the theoretical 



5'^This reduction of the \Fy^{s)\'^ peak value amounts to 13% O 30% for PQ = tttt, KK, and Kn, 
depending on the decay channel and the exchanged resonances accounted for. See also the analogue 
comparison for the decay width in Section 15.41 
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uncertainties. This resummation respects analyticity and unitarity only in a perturbative 
sense, a feature that can be improved [50] to hold to all orders as follows: 

1. The relevant phase-shift is obtained as the ratio between the imaginary and real parts 

of the vector form factor 6^^{s) = Im Fy^{s) /Re Fy^{s) . Fy^{s) stands for 

a form factor which can be obtained from Eqs. flMl) . fl26l) and fl27|) by taking out 
the exponentials and placing back the real part of the relevant loop functions in the 
denominator (see Refs. p9] and [52] for details). 

2. A three-times subtracted dispersion relation is used to resum FSI effects and the final 
form factor reads 

F-%) = exp !„,, + a.,s^ + - [^^^ d.'j^j^^,-^ | , (52) 

with Scut ~ 4 GeV^ [50 | I124[ 1125] and Sthr = ijnp + mQY. The subtraction constants 
«! and ^2 are to be fitted to data. 

This procedure will be followed in a future update of the program. Technically the 
method requires calculation of Cauchy Principal Value integrations which can be rather 
time-consuming and necessitates numerical stability checks. 

For the moment we consider FSI effects as they are given in Chapter 12. 4[ Eqs. (!23|) - (!27|) . 
To run the code with the FSI effects one has to fix FFVEC = 1 in value_parameter .f . 
Effects of FSI are presented in Table [6], last four lines are provided for comparison, and in 
Fig. [TOl they change by 14% — 32% the decay width, depending on the channel. Further 
results are collected in our Web page [28] . 

In the three-meson modes these interactions are neglected at the moment. We plan to 
include them in the future, at least for the r — )■ ttttttz/t- decays. 
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Comparison of Mass(1) of K- piO in chiannel tau- => K- piO nu_tau 



xlO-" 




Figure 10: Normalized events distribution vs. invariant mass of 7r^K~ pair: no FSI case 
corresponds to the red (darker grey) line, for FSI the green (lighter grey) line was used. 
Adjustment of numerical parameters is taken into account (Shift of K* mass, etc.). Ratio 
of the two histograms is given by the black one (only in this case the left side scale should 
be used). We have used MC-TESTER, Ref. [57], for preparation of the plot. Some of au- 
tomatically generated markings are explicitly left on the plot in this case. The program, 
MC-TESTER, is used also for plots of our Web page [28j . 
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Parameter 


Var. name 


Default 


[suggested range] 


Used in channel 


M, 


mro 


0.77554 


[0.770,0.777] 


1 


Mp 


mro 


0.775 


[0.770 


0.777] 


4-9,* 


Ma, 


mmal 


1.12 


[1.00 


1.24] 


5-9,* 


M^ 


mrhol 


1.453 


[1.44 


1.48] 


1 


M^ 


mrhol 


1.465 


[1.44 


1.48] 


4,5,6,* 


r.' 


grhol 


0.50155 


[0.32 


0.39] 


1 


Tp' 


grhol 


0.4 


[0.32 


0.39] 


4,5,6,* 


M,„ 


mrho2 


1.8105 


[1.68 


1.78] 


1,4 


r,« 


grho2 


0.4178 


[0.08 


0.20] 


1,4 


7 


coef_ga 


0.14199 


[0.077 


0.099] 


1,4 


5 


coef_de 


-0.12623 


[-0.035 


-0.012] 


1,4 


01 


phi_l 


-0.17377 


[0.5 


0.7] 


1,4 


4>2 


phi_2 


0.27632 


[0.5 


1.1] 


1,4 


Mk^± 


mksp 


0.89166 


[0.891 


0.892] 


2,3,7-9,* 


Mk*o 


mksO 


0.8961 


[0.895 


0.897] 


2,3,7-9,* 


Mk^ 


mkst 


0.8953 


[0.8951 


0.8955] 


2,3 


Mk^ 


mkst 


iMK*± + Mk,o) /2 






7-9,* 


rriK' 


mkst 


0.94341 


[0.9427 


0.9442] 


2t,3t 


Tk* 


gainma_kst 


0.0475 


[0.047 


0.048] 


2,3 


Ik* 


gainma_kst 


0.06672 


[0.0655 


0.0677] 


2t,3t 


Tk*, 


gainma_kstpr 


0.206 


[0.155 


0.255] 


2,3 


Ik*' 


gainma_kstpr 


0.240 


[0.120 


0.380] 


2t,3t 


Mk*' 


mkstpr 


1.307 


[1.270 


1.350] 


2,3 


niK*' 


mkstpr 


1.374 


[1.330 


1.450] 


2t,3t 


F 


fpi_rpt 


0.0924 


[0.0920 


0.0924] 


all,* 


Fk 


f k rpt 


1.198F 


[0.94i^ 


; 1.2F] 


3,4 


Fv 


fv_rpt 


0.18 


[0.12 


0.24] 


5-9,* 


Gv 


gv-rpt 


F'/Fy 


[O.xxF'^/Fv 


l.xxF'^/Fv] 


5-9,* 


Fa 


fa_rpt 


0.149 


[0.10 


0.20] 


5-9,* 


Pp 


beta rho 


-0.25 


[-0.36 


-0.18] 


5,6* 


IK-K 


gainma_rcht 


-0.043 


[-0.033 


-0.053] 


2,3 


IK-K 


gainma_rcht 


-0.039 


[-0.023 


-0.055] 


2t,3t 


Ov 


THETA 


35.26° 


[15° 


50°] 


7,8 



Table 4: Initialization parameters defined in file 

new-currents/RChL-currents/value_parameter .f : part 2, fit parameters. An educated 
guess for the variation of some of the resonance parameters is given. Energy units are powers of 
GeV. Channels identification numbers are defined in Table Ej Variables requiring rerun of pretab- 
ulation new-currents/RChL-currents/tabler/al/dalwid_tot_rhol_gauss.f are marked with 
*. The parameters coresponding to non-default currents of Kn modes (FFKPIVEC=0) are marked 
with ' . 
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Parameter 


Var. name 


Default 


[suggested range] 


Used in channel 


Ci - C2 + C5 


cl25 


0.0 




7-9 


25'4 + 95 


g4; g5 


-0.6 


[-0.64,-0.56] 


7,8 


92 


g2 


Mv 


[-33% 


+33%] 


7-9 


192-K^y/2Fv 


91-93 


gl3 


-2Mv 


[-33% 


+33%] 


7-9 


1927r2V2Fv- 


Ci - C2 - C5 + 2C6 


C1256 


3FvMv 
mn^V2F^ 


[-33% 


+33%] 


7-9 


d. 


d3 


M2 
647r2F2 


[-33% 


+33%] 


7-9 


di + 8(i2 — (is 


dl23 


0.05 


-50% 


+50%" 


7-9 


Ci + C2 - 8C3 - C5 


C1235 


0.0 


[-0.25 


+0.25] 


7-9 


C4 


c4 


-0.07 


free 


7-9 


£/4 


g4 


-0.72 


free 


7,8 



Table 5: Initialization part 3: odd- intrinsic parity sector. File 

new-currents/RChL-currents/value_parameter .f . Channels identification num- 
bers are defined in Table El The defaults, which follow Refs. [251 1121] . are needed to 
reproduce our figures. 
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No. 


Channel 


Width [GeV] 


Reference 


In new-currents/RChL-currents 

directory channel's current: 

file — )■ routine 



1. 


tt-ttO 


5.2441 


10-13 ±0.005% 


Subs. E 


1 


frho_pi.f ^ CURR_PIPIO 


2. 


'k'^K- 


8.5810 


10-15 ±0.005% 


Subs. E 


1 


fkpipl.f -^ CURR_KPIO 


3. 


n-R^ 


1.6512 


10-1^ ± 0.006% 


Subs. E 


1 


fkpipl.f -^ CURR_PIKO 


4. 


K-K^ 


2.0864 


10-15 ± 0.007% 


Subs. E 


1 


fkOk.f^ CURR_KKO 


5. 


Tl^Tl^Tl^ 


2.0800 


10-13 ±0.017% 


Subs. E 


D 


f3pi_rcht.f ^ F3PI_RCHr 


6. 


n%%- 


2.1256 


10-13 ±0.017% 


Subs. E 


D 


f3pi_rcht.f^ F3PI_RCHr 


7. 


K-n-K+ 


3.8460 


10-15 ± 0.024% 


Subs. E 


3 


fkkpi.f^ FKKPI* 


8. 


K^n-RO 


3.5917 


10-15 ± 0.024% 


Subs. E 


1 


fkkpi.f^ FKKPI* 


9. 


K-Tt^K' 


2.7711 -10-15 ±0.024% 


Subs. E 


1 


fkkOpiO.f^ FKKOPIO* 











*TheFi of formula ©. 


1. 
2. 
3. 

4. 


n-R^ 
K-K' 


4.0642 - 10-13 ± 0.005% 
7.4275 ■ 10-15 ± 0.005% 
1.4276 -10-1^ ±0.006% 
1.2201- 10-15 ±0.007% 


Subs. |2H 
Subs. |2.4| 
Subs. |2.4| 
Subs. |2j 


frho_pi.f -^ CURR_PIPIO** 
f kpipl . f -^ CURR_KPIO** 

f kpipl . f -^ CURR_PIKO** 
f kpipl . f -^ CURR_KKO** 



**FSI off 



Table 6: Collection of numerical results to be obtained from the demo-standalone. Pos- 
sible future extensions going beyond the published paper will be documented in this table 
in the paper version included with the distribution tar-ball. References to subsections may 
be replaced in the future by references to the forthcoming papers. Last column includes 
references to routines of the currents code. This table is complementary to TableEl Results 
for the case when Final State Interactions (FSI) switched off with the help of FFVEC = 
in file new-currents/RChL-currents/value_parameter .f are also given, then FFKKVEC 
= 1 was used. 
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